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MODULI SPACE FOR GENERIC UNFOLDED DIFFERENTIAL 

LINEAR SYSTEMS 

JACQUES HURTUBISE AND CHRISTIANE ROUSSEAU 


Abstract. In this paper, we identify the moduli space for germs of generic 
unfoldings of nonresonant linear differential systems with an irregular singular¬ 
ity of Poincare rank k at the origin, under analytic equivalence. The modulus 
of a given family was determined in [9] : it comprises a formal part depending 
analytically on the parameters, and an analytic part given by unfoldings of 
the Stokes matrices. These unfoldings are given on “Douady-Sentenac” (DS) 
domains in the parameter space covering the generic values of the parame¬ 
ters corresponding to Fuchsian singular points. Here we identify exactly which 
moduli can be realized. A necessary condition on the analytic part, called com¬ 
patibility condition, is saying that the unfoldings define the same monodromy 
group (up to conjugacy) for the different presentations of the modulus on the 
intersections of DS domains. With the additional requirement that the corre¬ 
sponding cocycle is trivial and good limit behavior at some boundary points of 
the DS domains, this condition becomes sufficient. In particular we show that 
any modulus can be realized by a fc-parameter family of systems of rational 
linear differential equations over CP^ with fc-|-l, k-\-2 or k-\-3 singular points 
(with multiplicities). Under the generic condition of irreducibility, there are 
precisely k + 2 singular points which are Fuchsian as soon as simple. This in 
turn implies that any unfolding of an irregular singularity of Poincare rank k 
is analytically equivalent to a rational system of the form y' = • y, with 

A{x) polynomial of degree at most k and Pi{x) is the generic unfolding of the 
polynomial . 


1. Introduction 


The local classification in the complex domain of germs of systems of linear 
differential equations, with a pole at the origin, 


y = 




• y, 


exhibits a qualitative shift as one goes from A: = 0 to fc > 0. In both cases, one 
can first go to a normal form by a formal gauge transformation g{x), that is a 
power series in x. Let us suppose for simplicity that the system is nonresonant 
i.e. the leading term A(0) is diagonal, with eigenvalues which are distinct (for 
k = 0, one would ask that they be distinct modulo the integers). One can perform 
a formal normalization to have A{x) diagonal, and a polynomial of order k. If we 
then proceed to the analytic classification, one finds that for k = 0, the formal 
classification is the same as the analytic classification, in the absence of resonance. 
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For k > 0, the situation is very different. The formal gauge transformation does 
not in general converge, and one only has analytic solutions on 2k sectors around 
the origin, with constant matrices (Stokes matrices) relating the solutions as one 
goes from sector to sector. If we further assume that A(0) = diag(Ai,..., A„), and 
that we have permuted the coordinates of y and rotated x to e*®a; so that 

(1.1) Re(Ai) > • • • > Re(A„), 

then the Stokes matrices alternate between upper triangular and lower triangular 
as one goes from sector to sector. Once one has fixed the formal normal form, 
the Stokes matrices provide complete invariants. While these can be thought of as 
generalised monodromies (e.g.,[I7]), the passage from the irregular case {k > 0) to 
the regular case /c = 0 is not immediate, since the monodromies for k = 0 have 
no limit at the confluence. This passage however is a natural one to consider, in 
particular when unfolding a system with an irregular singularity. Doing so sheds 
new light on the meaning of the Stokes matrices, and this has been studied in 
particular in m, E. in. 0. 

Indeed, one has a deformation from one to the other. Let Pe{x),e G be the 
generic deformation of x^~^^ as a polynomial of degree k + 1: 

(1.2) p^(x) = x^^'^ + H-h eia; + eo, 


and then consider a deformed system 


(1.3) 


A((i,x) 

Pe(x) 


For a generic value of e, the singularities are simple poles, and the classification, 
for a fixed formal form, is essentially the monodromy representation; at e = 0, 
and more generally, along the discriminant divisor A(e) = 0, one has higher order 
singularities and so the Stokes factors. 

In |14| . [S], the problem of studying the family, i.e., the unfolding of the origi¬ 
nal system was addressed. It involved the seemingly simple step of rewriting the 
equation as a coupled system 


(1.4) 

(1.5) 


y = A(e,x) ■ y 

X = pe(x) 


of a vector equation (1.41 and a scalar equation (1.51 in an extra variable t. 


One begins by analyzing the scalar equation, appealing to some quite elegant 
work of Douady and Sentenac [B]. One considers real trajectories (Im(t) = constant) 
of the scalar equation. Away from a real codimension one bifurcation locus, the 
results of jB] partition the e-space into Ck Douady-Sentenac domains, or DS domains 
Ss, all adherent to 0, where 


( 1 . 6 ) 


Ck = 


ill 

fc -I-1 


is the k-th Catalan number. The Ss can be extended to wider DS domains Ss, 
which retract to Ss , the union of which covers all values of e for which the singular 
points are all of multiplicity one, that is, the complement of the discriminantal 
locus. Each of these domains is contractible. In [B], these domains are referred to 
as sectoral domains. 
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Figure 1. The domains and associated Stokes matrices. 
In fact, for e 7 ^ 0 the domains are generically limited by spiraling 
curves in the neighborhood of the singularities: we have repre¬ 
sented the boundaries by segments for simplicity. 


On each Sg, and for each e € Sg, the x-plane is decomposed into 2k generalized 
sectors (see Figure fllb)), each adherent to two singular points of the scalar 
equation as in Figure [^bjl This generalizes the natural sectors of normalization for 
e = 0 ; indeed, the boundaries of the sectors, instead of all terminating at a point, 
terminate at different singular points, which are the various zeroes of (x). These 
zeroes are vertices of a natural tree, so that the singular point at e = 0 has in some 
sense expanded into a tree. 

Turning now to the vector equation, one has first (see M) a straightforward 
extension of the formal normal form to the deformed equation. Indeed, we recall 
that a linear system = A{x) ■ y, y G C", with irregular singular point of 

Poincare rank k, and leading order term with distinct eigenvalues, has a diagonal 
normal form 

(1.7) y'= ^^^(Ao-k AixH- \-Akx’^)-y, y&C^, 


where the Kj are diagonal matrices containing the (fc-|- l)n formal invariants of the 
system, and Aq has distinct eigenvalues. This extends to the deformed equation 
( 0 ): 


( 1 . 8 ) 


^^-^(Ao(e) -k Ai(e)a;-k-kAfe(e)a;'^') • y. 


J/eC”. 


For systems with a fixed formal normal form, one then wants an analytic classifica¬ 
tion: for this, one first fixes a DS domain Sg. For e € Sg, on each sector in the 
cc-plane, one has geometrically defined bases of solutions, defined up to an action 
of the diagonal matrices, which are such that in going from one sector to the next, 
the change of basis matrices generalize the Stokes matrices, and indeed exhibit the 
same alternation between upper and lower triangular. When one goes from Sg to 
Sg>, the decomposition bifurcates, and one obtains different matrices; we will see 
that these are related to each other by algebraic relations. 
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In [9], it was shown that these matrices (together with the formal invariants of 
( |1.8[ )) provide a complete invariant for the system. That is to say, if one fixes the 
formal normal form, the generalized Stokes matrices defined for each DS domain 
Ss determine the system: two unfoldings with the same invariants are analytically 
equivalent. The current paper addresses the realization problem: fixing formal 
invariants given by /c + 1 diagonal matrices Ao(e),..., Afc(e) depending analytically 
of e, and given Ck sets of 2k unfolded Stokes matrices, one for each DS domain 
Ss, depending analytically on e G Sg with the same limit for e —>■ 0 under which 
conditions does one have a system (1.41 corresponding to them? The answer turns 


out to be that the formal invariants and Stokes matrices define the same monodromy 
representations, and exhibit some natural regularity near the discriminantal locus. 
This is achieved in particular by realizing the equivalence class as the germ of a 
system defined over all of CP^, and exploiting compactness. In the course of doing 
this, we also discuss the development of normal forms, i.e. canonical representatives 
of equivalence classes. 

To prove the realization we proceed in two steps. The first step is to realize 
the modulus over each DS domain Sg in parameter space. We find that any formal 
data and Stokes matrices depending analytically on the parameters can be realized. 
When the singular points are simple, we obtain a Fuchsian system on CP^ with 
singular points at xi,..., x^+i, R, oo, where xi,..., Xfc+i are the zeroes of Pe, and 
R, oo are two fixed auxiliary points independent of e: 


(1,9) 






_ A{e) 

X — xe X — R 


This realized family is unique up to gauge transformations which are constant in x. 
As a connection, it is indecomposable; by, for example, diagonalizing A(0, 0) and 
normalising certain terms, we can make it unique. The residue matrix at infinity 
of ([T^ is simply A^{e) = - - ^(e)- 


The next step is to realize the modulus over a full neighbourhood of the origin 
in parameter space, which we can take as a polydisk Dp. For this, an additional 
condition is needed. Indeed, since the modulus characterizes families of systems 
of linear differential equations up to analytical equivalence, it is obviously a nec¬ 
essary condition that the realized families over the different DS domains be an¬ 
alytically equivalent over the intersection of DS domains. A necessary condition 
ensuring the equivalence is that the monodromy representations associated to the 
realizations over the different DS domains be the same (up to conjugacy as in the 
Riemann-Hilbert problem). The matrices conjugating the representations from one 
DS domain to the next must in addition form a trivial cocycle. This condition 
turns out to be sufficient, and can be expressed in terms of the modulus, i.e. the 
formal invariants and the Stokes matrices over each DS domain. Now, over each 


DS domain, we have realized unique normalized families of the form (1.91. Because 


of the normalization, they coincide on the intersection of the DS domains. Hence, 
we have realized the modulus over a polydisk Dp in parameter space minus the 
discriminantal locus A = 0. To extend the realization to the whole polydisk Dp, 
we first extend it to the regular points of A = 0 (where only two singular points 
coalesce). We then fill in for the remaining codimension two set of values of e using 
Hartogs’ Theorem. 
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The particular case where the system of Stokes matrices at e = 0 is irreducible 
is worth noticing: indeed, we can realize the data in a system 


( 1 . 10 ) 




■ y, 


which is Fuchsian when the singular points are distinct. 

Our realisation gives us a (local) normalisation. We close the paper with a 
discussion of normal forms. 


2. Preliminaries 

2.1. The scalar equation and DS domains. We recall in this section the results 
of the unpublished work of Douady and Sentenac [ 6 ] . As discussed in |9] the con¬ 
struction of the modulus in [S] was governed by the dynamics x(t) of the polynomial 
vector field 

dr f) 

( 2 . 1 ) x= — =p^{x) — , p^{x) = x''+'^ + €k-ix''~^ --keix-keo, 

at ox 

on CP^. It will suffice for the moment to limit ourselves to the set of generic values 

(2.2) So = {e : A(e) ^ 0}, 
where A(e) is the discriminant ofp^{x). 

The analysis centres on understanding the real flow lines in the x-plane, i.e. 
those given as the images of Im(t) = constant. Thus, we are restricting complex 
flow lines to a foliation of real lines; as such, the dynamics near the singular points 
Pe{x) = 0 has certain special properties, not shared with generic real vector fields on 
the plane. On Sq, each singular point xe has an associated eigenvalue =p'^{x(). 
Then 

• The point is a radial node if Af G K. It is attracting (resp. repelling) if 
A^ < 0 (resp. A^ > 0 ). 

• The point Xf, is a center if \i G iK. 

• The point xi is a focus if A^ ^ K U fM. It is attracting (resp. repelling) if 
Re(A^) < 0 (resp. Re(A^) > 0). 

Moreover such a vector field never has a limit cycle. 

Remark 2 . 1 . We emphasize that whether a singular point is of a-type (repelling) 
or uj-type (attracting) depends importantly on the family of real flow lines one is 
considering in the complex plane, in particular their asymptotic direction. On the 
complex line, there is no such concept. A linearized example suffices to illustrate. 
Indeed, if one is considering x = ax, one has the solution x = exp(at) over the 
complex line. Substituting t = , for 9j real and varying r through the positive 

reals, one has one solution exp(ae*®^r) which, when r —>■ -|-oo, spirals into the origin 
when 9i is chosen so that Re(ae®®'^) is negative, and one solution exp(ae*®^r) which 
spirals outward when Re(ae®®^) is positive and r —>■ -too. 

To understand the global structure of the real flow lines, the point x = oo serves 
as an organizing centre; indeed, the vector field Ve{x) has a pole of order fc — 1 there, 
and the system is structurally stable in the neighborhood of oo as e varies. 

Among the solutions Im(f) = constant, we have 2k separatrices at x = oo, alter¬ 
nately attracting and repelling (see Figure]^. On Sq, the dynamics is completely 
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Figure 2. The separatrices of the pole at oo. 


determined by the separatrices. Following the separatrices in from infinity, either 
backwards or forwards, one has: 

• For generic values of e one lands at repelling (Re(t) —)■ —oo) or attracting 
(Re(t) —)■ oo) singular points xg of focus or radial node type. For such an e, 
there exists trajectories joining the singular points two by two and choosing 
one trajectory for each pair yields a tree graph (see Figure [^(b), which we 
call the Douady-Sentenac tree. We denote by Ei the set of generic values of 
e in So for which there are no homoclinic connections between separatrices 
of oo. 

• The sets of generic e are separated by the closures of bifurcation sets of real 
codimension 1, where a homoclinic connection occurs between an attracting 
separatrix and a repelling separatrix of infinity: there is then a real integral 
curve flowing out from infinity in the x-plane and flowing back to infnity in 
finite time. On these bifurcation sets, the singular points can be split into 
two sets 1 1 and I 2 and 



(2.3) 


One sees this by integrating the form dt along a homoclinic orbit, and 
evaluating residues. When Ii is a singleton, the corresponding singular 
point is a center. 

The closures of these codimension one sets partition their complement in e-space 
into a certain number of connected components, which we call Douady-Sentenac 
domains, or DS domains. As one moves through the closure of the real codimension 
1 homoclinic locus, the limit points of attachment of the separatrices will change. 
These attachments are constrained, and the results of Douady and Sentenac tell us 
that there are 


ways of doing it (see Figure [^b)), thus dividing the complement of the closure of 
the homoclinic set into Ck open sets C Sq in parameter space. We will see how 
this happens below. 
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2.2. The sectors in x-space over C. The sectors will be defined in several steps: 
we first define sectors g , which cover C minus a set of measure 0, and then 
their enlargements > which cover C minus the zeros of Pe{x). We modify 

them later so that they are adequate on a disk. In this section we limit ourselves 
to sectors on C. 

Depending on the meaning one gives to the word “explicit”, the flow lines of the 
scalar equation are explicitly solvable. Indeed, one has a t{x), which globally is 
multi-valued with multi-valued inverse, defined by 

To fix ideas, when fc = 1, one can solve further: 


(2.5) 


t{x) 




log 


x+-^—eo ’ 


eo — 0, 

eo 7^ 0. 


For eg 0, the image of a disk D is the complement of a line of periodic holes 
located apart. At the limit when eg = 0, all holes but one disappear to 

infinity. 

More generally, since cx) is a pole of order k—1 for v^{x) (hence a regular point if 
fc = 1), it can be reached in finite time. Hence, the image of oo in cc-space consists of 
finite point(s) in t-space. Also, the time t is ramified at c» for fc > 1, since t ^ 
near oo. For e = 0, the image in t-space of a disk D in x-space is the outside of 
a disk on a fc-sheeted Riemann surface. For e 7 ^ 0, the map t is multivalued with 
the different images periodically spaced. The image of a disk is the complement 
of a countable number of periodically spaced holes placed on a branched fc-sheeted 
Riemann surface. The periods between the holes tend to 00 when e —)■ 0. (More 
details in jH] and below.) 

The fact that the integral curves are given in terms of a simple integral gives 
quite a lot of control over the behaviour of solutions. For e in a DS domain Sg in 
parameter space, we will get a decomposition of the complement in the cc-plane of 
the separatrices emerging from infinity in the complex cc-plane into 2 fc (generalized) 
sectors g , ordered cyclically at infinity in a clockwise direction (around infin¬ 
ity!) ^2,)L,Ss’ ^ 2 ,e,Sa> ^fe,e,Ss> as in Figure j^b). The clockwise 

direction at infinity will become the standard anti-clockwise direction when viewed 
from the finite regions of the plane. In a neighbourhood of a; = 00 the sectors are 
well defined: they are simply the complement of the separatrices in the bifurca¬ 
tion diagram at infinity. The sign ± reflects whether the sector has its outgoing 
separatrix on its right (-I-) or left (—) hand side, when looking outward from infinity. 

The question is then what happens as one extends inwards. Following the sep¬ 
aratrices inwards, each sector will be adherent to two vertices Xa, x^j, with x^j 
attracting and Xa repelling, and both points being singular points of Vc(x). The 
boundary of each sector will consist of three sides: 

• an attracting separatrix of 00 (i.e. going to infinity) emerging from a singu¬ 
lar point Xa of a-type (a source) of v^. in f-coordinate it is represented by 
a horizontal half-line (real values of time) starting from infinity to the left 
(corresponding to the singular point Xa) to a finite image in t of x = 00 ; 
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• a repelling separatrix of infinity going to a singular point x^j of oj-type (a 
sink) of some : in t-coordinate, it is represented by a horizontal half-line 
starting at a finite image of a: = oo and directed to the right; 

• and a curve (not uniquely defined) from Xa to x^j, corresponding to a real 
trajectory of Ve from Xa to x^. 

Thus the sector is topologically a triangle, though of course most of the time 
the boundaries of the sectors in a:-space will be spiraling curves. While giving 
the formulae for the boundary is impossible, because of the differential equation 
^ = Peix), the sector in t-space is fairly simple, and is given by a horizontal strip 
(see Figure [^. 



(a) Sector on C (b) Sector on a disk 


Figure 3. A strip in t-space whose image is a sector ^ 5 in cc-space. 


Construction of a pair of sectors in x-space. We examine the 

construction in a bit more detail, giving now a pair of sectors as the image of a 
(wider) horizontal strip in t-space. Let us choose a ^ for some ji, and consider 
the two separatrices on its boundary. These can be thought of as the image of 
one line Im(t) = constant in t-space, the union of two half-lines, with the first 
emerging at Re(f) = —00 from some singular point Xa, and arriving at some time 
to at infinity in the x plane; the second half-line starts at tg and then flows off to 
an x^ at Re(t) = -foo. Moving inwards in the x-plane near infinity in the sector 
j (and so downwards in the t-plane), the lines Im(t) = constant now become 
flow lines from Xa to x^)- This patterns continues downward in the t-plane, until 
one hits a line Im(t) = constant which again goes through a; = 00 at a time ti, 
which this time will be the boundary of a domain now viewed as a part of 

the plane above the <i-line. One has then created a horizontal strip in the t plane, 
bounded above by the line through and below by the line through ti, which 
corresponds to the union of two sectors 0“ we choose the demarcation line 
between the two sectors, somewhat arbitrarily, to be the image of the real line in 
f-space through {ti -|-to)/2 (see Figure]^. 

The above construction defines the sectors; it also does a bit more. Indeed, 
first of all, it pairs the sectors, with each sector g (and its enlargement g ) 



(a) Sectors on C (b) Sectors on a disk 

Figure 4. A strip in f-space whose image is the union of sectors 
s e s x-space. 
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Figure 5. The separating graph formed by the separatrices land¬ 
ing at the singular points and the curves 7 ^ (in dotted lines) used 
to calculate the r^. 


with a positive sign paired with one sector ^ 5 (and its enlargement e s ) 
with a negative sign. It also gives us, for each pair, a path joining the positive and 
negative sectors, simply as the image 7 ^- in a;-space of the segment joining to and 
ti; the curve 7 ^ is then a path out from infinity, going back to infinity (Figure]^. 
As Douady and Sentenac remark, the various 7 ^ do not intersect; their complement 
in a sufficiently large finite disk is a union of disjoint open sets, each containing a 
single zero of . Homotopically, our flow lines in a;-space give a diagram consisting 
of a disk with points representing the sectors g , 17)”^ g , 17^^ g , 17^^ g , , 

^ke s ’ that order, on its boundary, and non-intersecting chords (the 7 ^) joining 
the sectors which are paired. Dually, one has the Douady-Sentenac tree joining 
the various zeroes of . Douady and Sentenac show that this (“Douady-Sentenac”) 
invariant characterizes the DS domains, and it is the count of the possible diagrams 
that gives us a Catalan number. 

One also has more; the difference 



gives us a complex parameter in the upper half plane, and Douady and Sentenac 
show that the set of such parameters as one ranges through the paired sectors 
parametrizes the DS domain: 

Theorem 2.2. |6] There exists a biholomorphism : Sg ^ where El is the 
upper half-plane. The map is defined as follows: $s(e) = (D) ■ ■ • where 

f dx 
Jj.Peix)' 
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In particular, Sg is contractible, and so, simply connected. The set 

which we call the spine of Sg, corresponds to polynomial vector fields with only real 

eigenvalues at each singular point. 


Note that each Douady-Sentenac parameter Tj is computable as a sum of residues 
Tj — the zeroes xe surrounded by 7 ^. The fact that the 7 ^ 

partition the cc-plane into regions each containing only one zero of Pt{x) then tells 
us that each residue in turn is computable from the Tf. there is a bijective linear 
map between the Tj and the (satisfying = 0). The map 

between the space Sq of polynomials Pe with distinct zeroes and the quantities 
Pg{xj))~^ is many to one, a branched covering, but of course bijective to its image 
when restricted to a 5's. 

The radial rescaling on C and its effect on polynomials has some interesting 
consequences. The bifurcation diagram of the ODE (2.11 rescales naturally when 
changing {x,t) 1 —>■ {rx,r~^t), with r G K^, which induces the parameter change 

(2.6) (efc_i,...,ei,eo) (r^e^-i,..., r'^'+^eo) 


This suggests a natural “norm” in e-space given by 

(2.7) II e ||:= max (||eo|^^^''’^^^,..., |efc_i 



which is essentially the maximum norm of the roots of p^, and the following equiv¬ 
alence relation on the set of e: 


(2.8) e = (eo,..., €k-i) ~ e' = (ej,,..., e),_i) 3r e M+ : e'- = ^ej. 

This yields a conic structure in the parameter space, which is foliated by curves 
(equivalence classes) {(r^efe_i,..., r^ei, r*+^eo) : r S K'*'}, with fixed e. Both the 
sets Sg and the strata of the bifurcation diagram are unions of equivalence classes, 
which can be described by their intersection with the sphere || e ||= p for some 
positive p. 


Lemma 2.3. Let Sg be a DS domain. Then there exists C > Q such that, for each 
eG Sg with II e ||< 1, if^g{e) = (Ti(e),..., rfe(e)), then 


(2.9) 


lD(e)l > 


C 

RF' 


Proof. We consider the intersection of the subspace || e ||= 1 in parameter space 
with Sg. The quantity C = min||£||=i{|Ti(e)|,..., |Tfe(e)|} is bounded below on this 
intersection by a positive constant C. Indeed, we have seen that we have a bijective 
linear map relating the Tj to the allowed values of the quantities vi = p'{xi)~^ = 
— Xk)~^ (note that = 0 ), so that a bound for one is a bound for 

the other. On the other hand, the product ~ Xk)~^ is bounded below, 

as the roots are bounded above on the set we are considering. In particular, the 
product goes to infinity as one approaches Aq. 

Then, the rescalings of the roots of p^ and (2.6) yield |rj(e)| = ? from 

which the conclusion follows. □ 


Remark 2.4. The perceptive reader will have noticed a problem with our con¬ 
struction. When dealing with sectors in C, the construction has no limit when 
approaching the boundary of a DS domain Sg where t^ — ti G M. Also, the image of 
a disk Dfi in x — space is, for small e, approximately the exterior of disks of radius 
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1 

kR^ 


Even if to — ti £ H, where El is the upper half-plane, it may not be possible 
to pass a horizontal strip in the case of sectors on a disk (Figure^b)). Of course, 


since to — ti = 0(|j 


\-k 


) (see Lemma 2.3 below), for any nonzero argument of 


to — ti, it suffices to reduce e to be able to pass a horizontal strip. But we want 
a uniform description on a poly disk {|| e ||= p}. The solution will be to slant the 
strips. We come back to this in Section f2.4l 


2.3. Sectors at e = 0. We pause for a second to remark that when e = 0, the 
separatrices partition the Riemann sphere into 2k isomorphic sectors, rather like 
the portions of an orange, with vertices at x = 0. The point is that the sectors, 
with our normalisation of the equations, can serve as Stokes sectors. The real flow 
lines within each sector flow in tangentially to the separatrices at the origin (taking 
a blow-up z = re*®). Meanwhile, in t-space, the strips become half-spaces, as e 
moves to zero. 

This is, in some sense, the purpose of the construction: we are getting, for general 
e, sets which generalise the Stokes sectors. Unfortunately, we cannot do it uniformly 
in e; each DS domain has its own extension. 


2.4. Enlarging the sectors and the DS domains. We have thus, following 
Douady-Sentenac, described the generic locus Ei of systems with no homoclinic 
loops through oo as the union of a certain number of DS domains Sg, each isomor¬ 
phic to H^, and a boundary consisting of (the closure of) the set for which there 
is a homoclinic orbit joining two separatrices. Our approach for understanding the 
deformation of an equation will be to first study the deformations over the domains 
Sg, and then to make identifications along their mutual boundary. As things are 
now, one would then need to have information continuous up to the boundary, and 
then know how to interpret it. It will be simpler for us to extend the DS domains 
Sg somewhat to wider domains Sg, on the complement Eq of the discriminant locus, 
so that they overlap, and we can simply glue over open sets. As usual, we are only 
concerned with e small. 

Likewise, for a fixed e we will want to expand our sectors g , so that together 
they constitute an open covering of the complement of infinity in the x-plane. This 
is indeed quite simple: one simply widens the strips by a fixed amount p in t-space, 
everywhere except for a 77 -neighbourhood of the points to, ti on the boundaries of the 
strips, where one instead widens the angle of the intersection of the 77 -neighbourhood 
with the strip, from tt to tt -|- 20, increasing the angle byO< 0 < 7 r/ 2 on each side; 
one can choose, for neatness’ sake, that p = sin( 0 ) 77 . Let us denote our widened 
sectors by g , removing the tilde. 

A final issue is that we have partitioned the complex plane; we would also like to 
ensure that this partition restricts well to a disk, which was our original situation. 

Let us now consider the DS domains. 

Before extending the domains, let us highlight what we need. 

• The first thing that we need is a partition (now a covering) of the x-space 
by 2k sectors, with on the intersection of neighbouring sectors, separatri¬ 
ces extending out in some reasonable way from infinity and lying in the 
intersection of our extended sectors. 

• The characteristic feature of the DS domains is the DS invariant, which joins 
a sector g to the sector ^ g , joining our sectors pairwise; this must 
of course be preserved in our extension; likewise the dual invariant, the DS 



12 


J. Hurtubise and C. Rousseau 


tree joining the singular points. Note that both these invariants imply a 
partition of the singular points into a- and w-type. These invariants should 
extend to the Sg- 

The third thing that is required is some control of the aspect of the ap¬ 
proach of our “real” lines as they go into the singular points in x-space, 
or alternately the asymptotic direction one chooses to go out to infinity 
in t-space. One of two key reasons for this has already been given in Re- 
A simple singular point Xf, ol x = Pe{x) such that p'^{xt) ^ M 


mark 2.1 


can be of a-type or of a;-type depending on how we approach it. On some 
of the boundary of Sg, Xf, is a center since Re(p'(a:£)) = 0. But we could 
force it to keep the same a-type or w-type as in Sg by approaching it along 
appropriate spirals inside a sector g . Again, one wants our singular 
points to retain their type over the extended domains. Another reason, is 
that we will also want to define a flag of solutions, with decay rates, of our 
vector equation y = A(x, e) ■ y. This will again depend on the asymptotic 
direction in t-space. What we will want here is that the family of direc¬ 
tions chosen not deviate too much from the real one. Hence the a-type or 
w-type is not an absolute property of the point, but a relative property de¬ 
pending on a sector g adherent to xi. The same can be done with the 


other boundary parts of Sg corresponding to a homoclinic loop surrounding 
several singular points. 

• One final requirement, of a more technical nature, appears in the defor¬ 
mation of our strips as one approaches the boundary of the Sg. Since we 
are in the end dealing with a local problem, on a disk D in a;-space, we 
would like our sectors to partition the disk D in a clean way. As we are 
considering deformations around e = 0, we note that a neighbourhood of 
infinity in x-space translates into (half-) disks Dq, Di around the points to, 
ti in our horizontal strips in t-space. We would like, to separate our two 
sectors defined by the strip, that at least one of the real lines from Xa to 
Xtj not intersect the half discs. On the other hand one need not stick to 
real lines, and in fact what we really want is that the lines defining the DS 
invariant and the dual DS-tree stay inside our disk D. 

• On the other hand, staying in a finite region of the t-plane, i.e away from 
the singular points in x-space, one has a holomorphic differential equation; 
this is a flat connection, and parallel transport along any two curves that 
are homotopic on the complement of the singularities gives the same result. 
In other words, one can deform quite freely as long as one stays away from 
the singularities, and does not wander out towards infinity. 


Now that we know what we want, let us see what we get when we go to the 
boundary of Sg. The parametrization of the Sg places, quite conveniently, the 
discriminantal locus out at infinity in H^, with the width of at least one strip in 
t-space tending to oo. In particular, for e —>■ 0, each sector for e = 0 corresponds 
to the image of a horizontal half-plane in t-space. We are, as we said, interested 
in deformations of an irregular singular point, i.e. in e small. From the scaling 
action considered above in ( 2 . 61 , or simply by considering the relation defining x, 
this means that we can restrict our attention to DS parameters Ti,...,rfc with 
I Til > C II e 11“^. This fact is of some use in helping our constructions go through. 
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(a) Sectors on C 


(b) Sectors on a disk 


Figure 6 . A slanted strip in t-space whose image is the union of 
sectors g and ^ s cc-space. 


On the other hand, in real codimension one, as one is moving to and through 
a homoclinic orbit, it is easy to see that the width Im(T), where t = to — ti, of 
a strip is tending to zero, so that in the limit one has a real line in the t-plane 
passing both through to and ti ; the homoclinic orbit is then the image in a;-space of 
the real segment joining to and ti. Moving through changes the order of the strips 
and bifurcates some of the limit points at t = ±oo to other singular points in the 
cc-plane. When the homoclinic orbit surrounds a unique singular point, then the 
singular point bifurcates from attracting to repelling, or vice versa. 

The domains are defined by where the separatrices, the flow lines ending or 
starting at a; = (X) end up in the cc-plane. The idea for extending the DS domain 
is simple; basically, away from the singular points, we are dealing with one large 
complex (and multi-valued) orbit. There is nothing sacrosanct about the lines 
Im(t) = constant; one could equally well, for example, consider lines Im(e®^t) = 
constant, or even piecewise linear curves of that type; these are, after all, restrictions 
of the complex flow lines of the same complex equation. We will see that the angle 
/3 is important, but that there is some flexibility, as long as the angle f3 is kept small. 
Thus, for example fixing to = 0, if one is moving the DS parameter ti up to the 
real positive axis, one can also move the family of lines, for example to 
= constant, for a suitable (3{ti) so that the strip moves with ti and preserves the 
limit points of the separatrices, as Im(to — ti) becomes zero, and even beyond (see 



But Figure is a special case. In general, there is no reason why Xa and Xc^ 
would turn of the same amount and in the same direction. Hence, once we have 
passed to and ti, we may have to bend the strip again so that it makes an angle Pa 
on the left (resp. Pi^ on the right) with the horizontal direction. One would also 
have to deform in between as several strips come together; a precise description 
of the strips in t-space seems a bit complicated, but we will see that this will not 
be needed. What is important is that all half-strips having the same Xa or x^j 
limit point have the same slope, thus yielding in cc-space spirals with same rate of 
spiraling. Also they must cover a full neighbourhood of the limit point. Let us give 
the constraint when the limit point is Xa- Let B — where J is the set of 

indices of Tj corresponding to strips with limit point Xa ■ Then H = ± ^ . Hence 

the sum of the width (in the direction of B) of the half strips having limit Xa must 
be more than B so as to cover a full neighbourhood of Xa- 

Let us consider how our space decomposes both in the t and x-parametrisations. 


before deformation of Sg into Ss (see Figure^: 
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Figure 7. A slanted strip in t-space whose image is the union of 
sectors g and e s x-space. The slope on the right is 
chosen so as to approach x^j- 



(a) a;-space 



Figure 8 . The different regions in x-space and t-space. The DS 
invariant is in grey and the DS tree in green. 


• One has points ti in t-space, mapping to oo in x-space. These are sur¬ 
rounded by (diffeormorphs of) disks Di, with Di minus a slit mapping to 
the intersection of a pair of consecutive sectors g with the set |x| > R. 
One can cover the region |x| > i? by a finite number of the Di, i = 1,... ,m, 
with the lines joining ti and ^^(i) generating the DS invariant in x-space. 
Let all the Di lie inside a vertical strip \Re{t)\ < C. 

• Taking our singular points Xq, in x-space lying inside the disk |x| < 
i?/2, we have consecutive horizontal half-strips in the region \Re{t)\ > 2C 
covering disjoint disks Ba, in x-space around the Xq-, Xi^, with Re(t) < 
—2C for an a-point and Re(t) > 2C for an cu-point. The spiralling in or out 
of the strips on these disks is important, as indeed is the particular aspect 
(slope in t-space) at which it is done: indeed, the in or out tells us whether 
the singular point is attractive or repelling, an important feature of the DS 
invariant, and the aspect of the lines influences the various growth rates of 
our solutions to the differential equation. 

• Between the two, we have families of horizontal lines in t-space joining in 
x-space the Be to B^, with many of them leaving the disk |x| > R. The 
whole family is not important, and indeed all that is necessary is the isotopy 
class of one of the lines, say laui in x-space for each horizontal strip, as this 
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is what defines the DS-tree. One has in addition the segments defining the 
DS-invariant joining the circle \x\ = R to itself, and avoiding the Bj. Again 
it is only the isotopy that is important. 

This is the picture on a DS domain Ss- We want to extend this picture to a 
larger set, obtained (as we are on the generic locus of distinct zeroes of p^) by 
moving the individual zeroes of around in small balls. As well, we would like to 
ensure that even on Ss, we can ensure that the segments defining the DS invariant 
and the dual DS-tree stay inside. It is clear that this is possible: 

• One leaves the region |x| > i? as is. 

• With an isotopy one pulls in the segments of the DS-invariant and the 
segments of the DS tree so that they lie in the ball |x| < R/2. Moving the 
zeroes of p^ around, one can isotope the DS-invariant and the DS tree to 
follow. 

• One moves the balls Bi around with their zeroes, and modifies the slopes 
of the strips in t-space that map into them so as to satisfy our requirement 
of the points staying attractive or repulsive. 


This procedure highlights the need to define the new DS domains Ss as covering 
open sets in the etale sense; that is, as spaces above the complement Eg of A = 0 
(Eg has real codimension 2), which in fact contain the extra information of the 
attachment of the singular points to infinity and which retract to Ss- This is 
important when one branches around the locus A = 0: indeed, while the locus 
of the zeroes of p^ is the same, the pattern of attachment of the separatrices to 
the singularities can change; likewise, the segments of the DS-invariant or the dual 
DS-tree will follow the zeroes as they are isotoped. This was brought to the fore 
in the treatment of the k = 1 case in El; one sees there that the projection of 
Ss to the neighbourhood of A = 0 self-intersects, but with different attachment to 
infinity (see Figures and 13 later in the paper). This persists here: for generic 
points of A = 0 where exactly two points of = 0 coincide, a single Ss can be 
extended in order to cover the neighbourhood of that point in a ramified way. In 
addition, of course, neighbourhoods of points of A = 0 inside Eg can be covered by 
several Ss- 


2.5. Some conventions. We will be considering our linear o.d.e. problem from 
a slightly more abstract point of view, as lying on bundles equipped with singular 
meromorphic connections. As this involves changes of basis and so on, we pause 
to fix some conventions, some of which do not seem to be entirely uniform in the 
literature. 

A fundamental matrix solution to the o.d-e-, or of covariant constant sections- 
This will be a matrix X, everywhere non degenerate, with 

{d- A)X = 0. 

Then the coordinates of a solution vector, with respect to a basis of solutions is a 
constant column vector v, with 

X-v, 

representing the solution vector in our trivialisation. Likewise, two fundamental 
matrix solutions are related by a constant invertible matrix C 


Xi — X 2 ■ C 
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Changing gauge. The fundamental matrix solution gets multiplied by g on the 
left, so that X = gX-, simultaneously, the trivialisation of the bundle gets acted 
on by g~^ from the right. The matrix is still thought of as giving the same set of 
solutions, but given in a different gauge: we have 

dX-{dg-g-^+gAg-^)X = Q 


In particular, taking g = X~^ gauges the connection to zero. 

The Stokes matrices. Recall that the real flow lines at e = 0 partitioned the 
Riemann sphere into 2k sectors, all meeting at 0 and infinity. Around the origin 
or the boundary of a disk the unfolded sectors are ordered counterclockwise, for 


arbitrary DS domains and e, as 






l,e.S, 


Q+ 


n 


2,e,Ss 


n 


k,e,S, 


Note 


that this is the clockwise order if one turns around infinity. 

The theory of irregular o.d.e. gives us on these sectors (limited to a disk), 
fundamental matrix solutions X^ on asymptotic to a diagonal formal solution; 
on intersections, they are related by Stokes matrices: 








with a cyclicity convention in the indices that 0 = k. 
changes of coordinates in flat trivializations: 


Dually, they can be seen as 


^ 3-1 = ^3 ^3 


= C^^3 


(the v’s are all the same vector, but in different coordinates). 

We will, on each DS domain, be extending these Stokes matrices to generalized 


Stokes matrices 


s ) along with some extra diagonal matrices 


Normalization of the Stokes matrices. Geometrically, the Stokes matrices are 
defined up to the action of a certain number of diagonal matrices. Traditionally 


in the literature, the Stokes matrices Cf 




pL 


are normalized so that their 


diagonal coefficients are 1. We rather choose a normalization so that the product of 
the inverses of the Stokes matrices in the appropriate order yields the monodromy 
matrix along the loop going around all the fc + 1-singular points. One way to achieve 
this is to take all Stokes matrices g^ with diagonal coefficients equal to 

1 , except Ss’’ which incorporates the formal monodromy exp (—27riAfc(e)), with 
—Afc(e) the residue matrix at infinity in (1.7). 

Monodromy. We follow the conventions of uni. The monodromy is defined by 
taking a fundamental matrix solution X, continuing it around a loop 7 , and relating 
the result X^ at the end to the beginning X = Xf, hy Xe = Xb ■ M , with M a 
constant matrix. On coordinates Vb (resp. vf) of a solution in a basis given by 
the columns of Xb (resp. X^), we have Vb = M ■ v^. This convention gives a 
representation of the fundamental group in the sense that 


-^71 *72 — -^72 ' -^ 71 ) 

where the composition 71 *72 of two loops is obtained by going through 71 first, then 
72 . Considering the monodromy around the origin, starting from Xjf on , the 
extension of XJf on is X^ (C^)”^, the extension of X^ on Tlf is Xf ..., 

giving the monodromy around the irregular singularity 

M = ■ - • {ctr^ ■ (cf )-i. 

This coincides with the conventions of m. More generally, if we have, moving from 
sector A to sector B, matrices Cab relating flat bases Xb = XaCab, then if along 
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a loop, if we go through sectors Ai,A 2 ,Ak, Ai in turn, the monodromy will be 
M = CAi,Ak ■ Ca^A^-i ■ ■ Ca2Ai- 


2.6. Flags of solutions and generalized Stokes matrices. Let us now turn 
our attention to the vector equation (1.4|: 

y = A{€,x{t)) -y. 


We will in particular be interested in the behaviour of solutions as x{t) goes to the 
limit points Xa,Xaj, i.e. as t goes to ±oo along our real lines. The dominant term of 
the equation is then A{e,Xa) or A(e,a;(^), which we denote j4(e,a:*). We then have 
the approximate equation: 

y = A{e,Xt) ■ y. 

As we are deforming from a generic irregular singularity, we can suppose 


A{e,x^) = diag(Ai(e),..., A„(e)), 


and that 


(2.10) ReAi(e) > Re A 2 (e) > • • • > ReA„(e). 

The approximate equation has solutions along each real line j/*(f) = (0,.., 0, exp(Ait), 0, ...0). 
The solution grows exponentially more slowly than j/® as t —>■ +oo, and grows 
exponentially faster as t —?► —oo. We thus have two flags of subspaces in the solution 
space, defined by growth rates: the first, given by behaviour as t —>■ +oo, defined 
by Wi{xaj) =< 2/"“®“'"^, .., y" >, and the second, given by behaviour as t ^ —oo, 
defined by Wi{xa) =< y^, ■•,y® >, with 

Wi(x,^) C W2(a::^) C ... C 


W„{xo,) D D ... D lFi(a:„). 


These two flags are transversal for sufficiently small e (see Lemma 2.6 below): 

Wi(x^) n Wn_i{Xa) = 0. 


At a generic singular point, for generic e, the actual behaviour of solutions is given 
precisely by the formulae above: the system admits a gauge transition to a normal 
form. Such a gauge transformation is not, however necessary, and indeed these flags 
are remarkably robust as one varies in e. One can see that this is reasonable in that 
the flag manifold is compact, and so there is at least a convergent subsequence as 
one moves around in e. In any case, the existence of the flag, and its continuity as 
one varies e (even moving out of Eg) is given by the following theorem of Levinson 
[TB] (see Coddington and Levinson |3], Theorem 8.1, p. 92; see also the discussion 
after Theorem 5.3 of 0). 


Theorem 2.5. Let a system of linear differential equations of the form 
(2.11) y = (^Ao(e) + A(e, t) + P(e, t)^ ■ y, 


he given on the real line, for which Ag is diagonal, with distinct real parts of the 
eigenvalues, A(e,t) is also diagonal, with limit zero at t = oo and 





dt < oo, 



( 2 . 12 ) 


\P{e, t)\dt < oo. 
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Then, setting A^(e, t), £ = 1,... ,k, to be the successive eigenvalues o/Ao(e)+A(e, t), 
there exist Iq S (0, cxd ) and solutions of the system for t S {Iq, oo) with 

rt 


lim ■ exp ( - / Xe{T)dT ) = ve{e), 

V Jto J 


for Vi(e) a non-zero eigenvector o/Ao(e) corresponding to A^(c»). //Ao(e), A(e,t) 
and P{e,t) depend continuously (resp. analytically) on e over compact sets in t- 
space, with the integrals in (2.121 uniformly bounded, then the solutions can be 
chosen depending continuously (resp. analytically) on e. 

In the cases of concern to us, the flags are varying continuously, and indeed 
analytically within the DS domains. In particular, we have 


Lemma 2.6. For e small within a DS domain, the flags Wi{xcj), Wn-iixcy,) are 
transverse in each strip. 


Proof. The flags are transverse at e = 0. Indeed, the standard results of the theory 
of irregular singularities (see for instance [TU]) give us on each sector a basis of 
solutions asymptotic to the standard basis of solution of the formal normal form; in 
these bases, the flags Wflx^), Wn-i{xa) are simply the standard flags, and indeed 
are transverse. □ 


Next, within a sector (i.e., a strip in t-space), the flags are constant, i.e. inde¬ 
pendent of the real flow line chosen. The reason is that one is solving an o.d.e. on 
the complex plane; one has a flat connection. Thus, to get the asymptotics of the 
flow lines on Im{t) = Ci, one can solve along Im(t) = cq, and then just integrate in 
the imaginary direction "out near infinity", on a segment of length |ci — co|, where 
the coefficients of the equation are essentially constant, and the equation preserves 
the flag. 

One thing to which the flags are sensitive, however, is the asymptotic direction 
of our flow lines in t-space. Indeed, if one just looks at the formal normal form, 
and takes the solutions y^{t) = (0,.., 0, exp(Ait), 0, ...0) for t = e*®s, s G K"*", it is the 
ordering of the real parts of A^e*® (assuming they are distinct ) which determines 
the flag. Unless 9 is suitably small, this is not the same as the ordering of the 
real parts of Xi. In particular, we must impose the constraint that the deformation 
angles be kept small in our construction of the extensions of the DS domains in 
the previous section. This is possible as long as we keep the parameters of our 
construction within a suitable range, which we can. 


Flags and bases; generalized Stokes matrices. On each sector, one has two transverse 
flags Wi{xuj), Wj{xa), which gives one dimensional intersections 

n Wn—i-\-l{^cf) —^ ej > . 

and so a basis Atr =< ei,e 2 ,...,e„ > defined up scale, i.e. to the action of the 
diagonal matrices. Now let us consider two consecutive sectors, sharing a separatrix; 
as such, they also share a flag. For instance, ,, share a separatrix adherent 

to a unique singular point Xa of a-type; along these we have a shared flag Wi{xa), 
and so their bases on ^ and aD^ on are related by an upper 

triangular matrix : 

Atr+ = Atr-_^ ■ 
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Similarly, and share a separatrix adherent to a unique singular point x^j 
of oj-type; along these we have shared flag Wj{xa;) and so their bases are related 
by a lower triangular matrix g : 

^tr- = Atr+ ■ 

Finally, consider and where a is the (order two) permutation represent¬ 

ing the Douady-Sentenac invariant; these share two singular points, one of a-type, 
one of w-type. The sectors and ^ then share both flags, and so the matrix 
e’ '^®l^ting the adapted bases, called gate matrix, will be diagonal: 

= Atr+ ■ 

These matrices will be our generalised Stokes matrices and gate matrices. 

2.7. Normalisations. We remark that there is freedom on each sector of an action 
by the diagonal matrices, as the intersections of the flags do not specify a scale. The 
Stokes matrices Cab associated to passing from sector A to sector B then get acted 
on by DaCabD^^. This is somewhat different from the classical Stokes matrices, 
which are usually defined as being strictly upper triangular, i.e with ones on the 
diagonal; this arises from the asymptotics of the equation. Here we allow an arbi¬ 
trary diagonal, but then quotient by the action of the diagonal matrices. We note 
then that our Stokes matrices contain the information of the formal monodromy 
around the singularity. 

Various normalisations are possible. One would like one such that the Stokes 
matrices g and g have limits when the confluence of singular points occur. 
We choose a normalisation that accomplishes this by taking all the diagonal parts 
of g^ and equal to the identity, except for the diagonal part g^ which 

is taken as exp(—27rzAfc(e)); this maintains the compatibility of our Stokes matrices 
with the formal monodromy. On the DS domains, one also needs that our Stokes 
matrices carry the monodromy of the formal normal form at each regular singular 
point; the product of the inverses of the diagonal parts of a sequence of the g 
or and of the corresponding to the sectors touching the singular 

point should be the formal monodromy. For example, if all the diagonal parts of 
the g ^, corresponding to the singular point are the identity, and there 

is only one gate matrix, then that gate matrix or its inverse must be the formal 
monodromy. In general, as the singular points are arranged in a tree, one sees that 
one may inductively move down the tree to compute the gate matrices in terms 
of the formal normal form. These normalised gate matrices C(^cr(j) e often have no 
limit when e —>■ 0, though the singularity disappears when we pass from this Stokes 
matrix picture A to picture B, as explained in Section [2.8| below. 

In what follows, we will be referring to normalised (generalised) Stokes matrices 
and Cl,^g^ ; we will suppose that the concomitant gate matrices have been 
computed from these and from the formal normal form. 

2.8. Building from the formal normal form; bundles on the disk or on the 
plane and Stokes matrices. We can consider our equations in a slightly more 
abstract form, as bundles equipped with singular (meromorphic) connections over 
the complex line. This allows us to choose different open sets covering the line, 
and to trivialize the bundle differently on these open sets. In our construction, we 
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will consider as basic open sets, for e G Ss, the sets g . There are three basic 
pictures we want to consider: 

A) Trivial connections d on each g , and constant transition matrices, defin¬ 
ing a bundle on a punctured domain (either C\ {cci,... Xk+i} or D\ {a:i,... 0 :^+ 1 }). 

Here we have on each sector o a bundle with the trivial connection, and 

Jit,Os 

transition functions which will be automorphisms of the trivial connection (i.e. 
constant matrices; these will be the Stokes matrices). The bundle on each sector is 
equipped with two constant flags Wi{xa) associated to the a- singularity adherent 
to g ^, and Wi(xi^) associated to the w- singularity adherent to 0+^ ^ 5 ^ ■ 

One has (covariant constant) trivializations compatible with the flags, in the sense 
that the i-th element of the basis lies in the intersection of Wi{Xa) and Wn-i+i{Xi,j). 

As one moves from sector to sector, the bundles are glued by constant transition 
functions, which are our generalized Stokes matrices. The intersections of two 
sectors, which we call intersection sectors, are of three types: 

• Sectors the intersection of containing a separatrix emerg¬ 

ing from Xa, a point of a-type; along these we have a trivialisation AtrJ_i 
on g and trivialisation Atr'^ on related by the upper triangular 
matrix g ^: 

Atr+ = Atr-_^ ■ 

• Sectors the intersection of n containing a separatrix converg¬ 
ing to Xuj, a point of w-type; along these we have trivializations related by 
a lower triangular matrix g : 

Atr~ = Atr+ ■ 

• Gate sectors gj the intersection n g, where tr is the per¬ 

mutation representing the Douady-Sentenac invariant; the gate sector is 
adherent to two singular points, one of a-type, one of w-type. The sectors 

and g intersecting in the gate sector share both flags, and the 

trivializations are related by a diagonal e ■ 

Atr-^^) = Atr+ ■ 


B) Diagonal connections = d— on each g , and transition matrices 

tending to the identity at the singular points, defining a bundle on the plane C. 
This picture is obtained from picture A by applying on each sector the gauge 
transformation given by the diagonal fundamental matrix of solutions to the 
formal diagonal normal form, (I. 8 I. (The fundamental solution given above is 
multivalued, and so one can choose different determinations on each sector.) The 
Tj*g preserve on each sector the flags Wi{Xa) and Wiix^S), which now however 
have geometric meaning, as they correspond to the decay rates of solutions along 
lines e^^t,t G K’*'. The transition matrices gauge transform to automorphisms of 


Along our intersection sectors: 


+ ^-l 


In H^g, we have matrices ~ 

of trivialisation from to : 


representing a change 
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These preserve the flag Wi(xo:), are upper triangular, with constant diago¬ 
nal terms and decaying off diagonal terms as one goes to the a-type singular 
point. 

In we have matrices 


(2.14) 


Btr^ = BtrJ 


(~iL 

'-'3,e,Ss 


These preserve the flag Wi{x^), are lower triangular, with constant diagonal 
terms and decaying off diagonal terms as one goes to the w-type singular 
point. 

> we set Ctuv. = 


In 


(2.15) 


= Btrt ■ cf 


it is natural to take the Fj^ as extensions of F^^ when 
± 


VO) - 

is diagonal. 

With our normalisation 
starting in g and crossing all g counterclockwise on a circle around the 
singular points. Then the gate matrices C^cr{j) e exactly compensate for the rami¬ 
fication of F^^ and = id. All our transition matrices have finite limits at 

the punctures, and so the bundle extends to the punctures. 

C) A globally defined singular connection d— with only one trivialization 

over a disk or C; the formal normal form of is . Again, associated to 

each limit point Xa or x^^ there will be (covariant constant) flags Wi(xa), Wi^x^j) 
on each sector defined by growth rates as one goes to the singular points along the 
paths Im(e*®t) = constant. In the x-plane, these paths are generically logarithmic 
spirals. On each sector, there is then a pair of (covariant constant) flags Wi{xa), 
Wi{xi^). The flags are transversal for sufficiently small e, by Lemma 2.6 This gives 
a unique (up to the action of diagonal matrices) basis of flat sections on each sector, 
with the i—th element of the basis living in Wi(xct) H Wn-i+i(xuj). We denote the 
fundamental matrix solution on each sector by g . 

Points of view A) and B) are easily seen to be equivalent; the point of this 
paper is to show that they are equivalent to C). As noted above, the gauge trans¬ 
formation from A) to B) is given by the F^^. Relating A) and C), one has the 

gauge transformation g . The gauge transformation from C) to B) is then 
fj± — z?± /' wi 




The monodromy of the connection on a path 7 around several singular points is 
given in different ways: in version C), one integrates the connection, as usual, and 
applies the convention above. In version A) it is given by the product of matrices 
g , g , and e their inverses, taken in the reverse order one meets 

the corresponding intersection sectors as one moves along 7 . In version B), it is a 
hybrid, a product in the right order of the matrices CW, g , g , ^ or their 

inverses and of the parallel transports by the diagonal connection on the intersection 
of 7 with each sector. 


3. The realization over a DS domain 


In this section we often drop the index Sg, as we will be working with a fixed 
DS domain. Also, we are concerned with e small, so in fact we only need to push 
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through our constructions on the intersection of Sg with a polydisk Dp around the 
origin in e-space. 

3.1. A bundle with connection at e = 0. For fixed e S Sg, the passage from 
versions A) or B) above to version C) is fairly well known, and the construction 
depends analytically on e G Sg. It amounts to finding the necessary gauge trans¬ 
formations on each sector to make the cocycles ^ trivial. Various 

techniques do this; we refer to HH]. This realizes the systems locally over a DS do¬ 
main. However, we will want to glue these realizations over DS domains to obtain 
a global family of systems for e in a neighbourhood Dp of 0. This glueing involves 
using the action of the gauge group to glue the systems. Gauge transformations 
over C form an infinite dimensional group; we would like to reduce the degrees of 
freedom somewhat, and we do this by compactifying. 

One would hope to realize the systems over Sg as singular Fuchsian connections 
on a trivial bundle over CP^, by adding in an extra singularity at infinity carrying 
the required monodromy. This is generically possible, but not always, as was shown 
by Bolibruch [3] and Rostov m- On the other hand, if we allow two singularities, 
then we can do it, at least for e in a small neighbourhood of the origin. 

We first realize the system for e = 0 as a system on a trivial bundle over CP^ 
with an irregular singularity at the origin, and two Fuchsian singularities, one at 
infinity, and one at some point at a large distance R from the origin, along the 
positive axis. This will reduce the gauge transformations to constant matrices in 
Gl{n,C). We would like the system to be rigid, in a suitable sense: 

Definition 3.1. A system of linear differential equations y' = A{x) ■ y is indecom¬ 
posable if it cannot be gauge transformed to a block diagonal form. 

Definition 3.2. A bundle with connection is reducible if it admits a nontrivial 
subbundle invariant under the connection. Otherwise, it is irreducible, i.e. the 
connection cannot be conjugated to a block triangular form. It is then also inde¬ 
composable. 

Remark 3.3. We show in Lemma 3.10[ below that any indecomposable connection 
can be normalized to a unique normal form. Then the same will follow for an 
irreducible bundle with connection. 

Remark 3.4. Choice of base point and trivialization. We will consider the 
monodromy of the connection at a base point Xf, = 3i?/4 , with some paths £R,ioo 
around the singularities at R, oo, whose product iuioo is homotopic to an anticlock¬ 
wise loop around the origin. We choose as trivialization of the bundle the canonical 
coordinate on flfg in which the connection (in picture A) is the normal form. We 
will choose changes of trivialization (to pictures B and C) which are the identity 
at the base point, so that the monodromy computed from this base point stays the 
same. 

With respect to a global trivialization (picture C), we will obtain a connection 
of the form 


(3.1) 



Let Aoo = —A)^ — Ar be the residue matrix at infinity (which vanishes when oo is 
a regular point). 
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Theorem 3.5. Suppose given 

• formal invariants given by diagonal matrices Ag,... Afc such that Ag has dis¬ 


tinct eigenvalues satisfying (1.11, determining a formal normal form (1.8 I, 
normalized invertible (Stokes) upper (resp. lower) triangular matrices , 
(resp. Cf ), j = 1,... ,k determining a monodromy 


(3.2) 


M = ■ (Cto) 


u \-l 


(Cton ■ (C^o) 


U \-l. 


• a generic matrix M^o representing a conjugacy class with distinct eigenval¬ 
ues, with all entries nonzero, close to the identity, and such that ■ M~^ 
has distinct eigenvalues. 

Then there exists a globally trivialized irreducible rational linear differential system 
(3.11 on CP^ with 


formal normal form (1.71 at the origin. 


Stokes matrices and Cj", j = 1,..., k, 

• monodromies in the global trivialization M{£ao) = Moo around 

in, loo, both with distinct eigenvalues. 

The automorphisms of the system are multiples of the identity. Furthermore, 
acting by an automorphism of the formal normal form (the constant diagonal ma¬ 
trices ), which transforms both the Stokes data and our rational system, it is possible 
to normalize the coefficients of n — 1 suitable monomials of entries of B{x) (or of 
n — 1 non diagonal entries of either the monodromies M{£ao)) to 1, thus 


leading to a unique normalization of (3.11 for each equivalence class under auto¬ 
morphisms of the Stokes data. 


Proof. Let us first build a bundle in picture B of Section |2.8| above, in the back¬ 
ground of the diagonal connection. We take our bundle on a neighbourhood of the 
origin, defined on sectors H^g, with transition matrices Cj^g,C^g. By the classi¬ 
cal result going back to Birkhoff [^1, this gets realized in picture C as a singular 
connection of Poincare rank fc on a disk D/j of radius R. 

From the base point Xb = 3i?/4 with the trivialization of the bundle fixed above, 
the monodromy M of the connection on the circle of radius 3i?/4 around the origin 
starting in is given by ( |3.2[ ). Now, build a holomorphic connection on the 
annulus Ar = {r G (i?/2, 3i?)}, with a singularity of Fuchsian type at x = i?, whose 
monodromy along the circle r = ^ (resp. r = 2R) is M (resp. the identity). (Note 
that it could happen that we start with x = R regular when M is the identity.) We 
glue this to our singular connection on the disk of radius R. The result is a bundle 
with a connection over the disk of radius 3i?, with an irregular singularity at the 
origin, and another Fuchsian singularity at x = i?. It has trivial monodromy around 
the boundary. Now glue this to the bundle on the disk 0.^0 = {r > 2i?} with the 
trivial connection. The result is a global bundle with connection on CP^ with two 
singularities, one at the origin, and one at x = i?. Any bundle on CP^ decomposes as 
a sum of line bundles; let this one have holomorphic type 0 (ai)©(P(a 2 )©...© 0 (a„). 

We now switch trivializations for a while, to the standard trivializations on 
®?=i^(oj)j®i Si CL 2 Si ■■■ S o-n- Let A be our connection in this trivialization, 
and let z = -. Let S = Im(x“’) = diag{z~°‘') be the transition from Uoo, a 
neighborhood of oo (resp. 0) in x-space (resp. z-space), to Uq = {x oo}. On 
Uoo, the connection is represented by a holomorphic matrix A(z) = (Ay(z)). A 
Schlesinger transformation on Uoo given by S' = Im(z““*) modifies the bundle to 
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the trivial one, modifies the connection matrix by SAS~^ + dS S~^, and changes 
the off-diagonal terms of the connection to SAS~^, i.e, Aij changes to AijZ°'^~°'% 
thus introducing poles potentially of order greater than 1 below the diagonal if Aij 
is suitably generic. 

We thus want to normalise the Aij by killing the terms of order less than ai — aj 
in z in the lower triangular terms Aij before conjugating by S. 

To do this we use the automorphisms g of (BiO{ai). These are also given by 
invertible lower triangular matrices, which in the C/qo trivialization are polynomials, 
with gij of degree Ui — Oj, representing a section of Hom(0(aj), 0(ai)). Indeed, 
because of the constraint on the degrees of the gj, then SgS~^ is again a lower 
triangular invertible polynomial matrix in x. 

Let us consider automorphisms g{z) such that g(0) = Id, i.e. g{z) = Id -f N{z), 
with N{z) = 0{z) lower triangular nilpotent. The automorphism changes the 
connection by 

A^B = gAg-^+g'g-^ 

One wants to get rid of terms in Bj of degree less than rriij = Oi — Oj for i > j; in 
other words, solving for some of the Taylor series of the differential equation above. 

Since N is nilpotent, then g-^{z) = Id - N(z) + (z) --f (-l)"-iAf"-i(z). 

Let B{z) = (Bjlz)), N{z) = {gij\z)) and E{z) = g{z)A{z)+g'{z) = (Eij^z)). And 
let Bj,g'^j,Aj,E^j be the terms of order fc in z of Bj^z), gij{z), Aij{z), Ej^z) 
respectively. Then Ej = Aj + {k + + QijiA,g), and so 

^ij = + {k + ^)9ij^^ + Qiji^^g)- 

where Qj and Qij are polynomials in the A^^ g, and g^i,, with j' < k. Hence we 
can solve successively degree by degree the equations corresponding to asking that 
the terms of Bj of degree m be equal to zero, starting from m = 0; solving for Bk 
“uses up” the form of the equations is such that all equations corresponding 

to putting some terms of degree m in some Bj with i > j are uncoupled, and hence 
can be solved independently for a fixed degree. 

Applying then the Schlesinger transformation S to B makes the bundle trivial 
and introduces a Fuchsian singularity (via the S~^dS term) at oo. 


Remark 3.6. We could have achieved the same purpose by using a polynomial 
gauge transformation on Uq and the permutation lemma of Bolibruch (Lemma 
16.36 of |10| ). so as to bring the singular point at infinity to be Fuchsian. 


Remark 3.7. Note that, generically, if we choose the right residues for our Fuch¬ 
sian singularity at R, the bundle is already trivial and we do not need to perform the 
Schlesinger transformation. In that particular case, the residue matrix at infinity 
Aao.o vanishes. 


We have built a bundle on CP^ with two sets of trivializations: the first, version 
B, has as open sets the sectors and the disk Doo = {\x\ > i?/2}; it has the transition 
matrices CWg,(7^Q and T between 11)1^ and the annulus A^. In this trivialization, 
the monodromy around R is M, and the monodromy around infinity is the identity. 
One also has a global trivialization, picture C, where our connection at e = 0 is 
of the form (3.1); we can choose this trivialization so that the leading term of the 
connection is Ag. So far, though, the monodromy around infinity is still trivial even 
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though infinity can be a singular point, and the bundle-connection pair might have 
non-trivial automorphisms. 

We would like 


• that our pair (bundle, connection) be irreducible, since then a further nor¬ 
malization will allow bringing it to a unique form (see Lemma 3.10 below); 

• that the point at oo have diagonalizable monodromy with distinct eigen¬ 
values; 

• that the singular point a,t x = R have diagonalizable monodromy with 
distinct eigenvalues. 

To do this, we deform the connection in picture B, keeping the same Stokes ma¬ 
trices, but modifying the monodromy around infinity from trivial to Mao, while 
modifying the monodromy around R in the opposite direction, so that the mon¬ 
odromy along the circle of radius R/2 stays constant. This uses Lemma 3.8 below. 

One can arrange that the monodromy around R is also of the desired form, by 
Lemma 13.81 below. 

Since Mao has all its coefficients nonzero, the system is irreducible (proof in 
Lemma 3.9 below). 

A key point is that a small deformation of a trivial bundle on CP^ remains trivial. 
Hence, for Mao close to the identity, our result is still a trivial bundle. Passing to 
our picture C, of a global trivialization, we get our desired connection. 

One can normalise the coefficients of the connection as in Lemma |3.10[ thus 
ending the proof of the theorem. □ 


Lemma 3.8. Any invertible matrix B can be written as a product B = C 1 C 2 of 
two invertible diagonalizable matrices Ci and C 2 with distinct eigenvalues. One of 
the Ci may be chosen in a Zariski open set. 

Proof. The map P : GL(n,C) —)■ GL{n,C) defined by P(Gi) = Cf^B is holomor- 
phic (even algebraic) and invertible. Let V C GL{n,C) be the subset of diago¬ 
nalizable matrices with distinct eigenvalues. It is the complement of an algebraic 
subset (where the discriminant of the characteristic polynomial vanishes). Then 
P{'D) n I? is the complement of an algebraic set in GL{n,C), hence nonvoid of full 
measure. □ 


Lemma 3.9. Suppose that the monodromy Mao has all its coefficients nonzero 
in the trivialization of the bundle over in which the connection is given by 
the formal normal form ([13 (a “picture B” trivialization). Then the system is 
irreducible. 

Proof. The point is then that from the point of view of the singularity at the 
origin, any decomposition must respect the geometric structures at the origin, for 
example the flags of growth rates. In the trivialization in question, these flags are 
the standard ones; this forces an invariant subbundle to be generated by a subset 
of the basis vectors. On the other hand, applying the change of basis to bring 
these basis vectors to be the first k vectors, which is a permutation matrix P, the 
transformed PMaoP~^ still has all its entries non-zero, and so does not leave the 
subbundle invariant, as it would have to be block triangular to do so. □ 

Lemma 3.10. We consider an indecomposable linear differential system (13, 
where Aq is diagonal with distinct eigenvalues. Then there exists n — 1 distinct 
pairs (ii,jt), with ie je, i = 1,... ,n — 1 and exponents € {1 ,..., A: -|- 1}, 
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such that the equation is conjugate by means of an invertible diagonal transfor¬ 
mation to a unique system with the coefficient of the monomial of 

normalized to 1. The only automorphisms of this normalised unique system are 
scalars cin for some c S C*. 

Proof. The normalization is done by action of invertible diagonal matrices. In 
view of the fact that the matrices c/„ are symmetries of the system, we limit 
ourselves to diagonal matrices D = diag{di,... ,dn) with di = 1. In order to 
be able to normalize n — 1 coefficients, we need to show that ii D ^ D' are two 
such matrices such that the conjugates of — R)y' = B{x) • y by D and D' 

are the same, then the system is decomposable. Let {!,..., n} = lU P, where 
di = d[ if and only if i G /. Then both I and I' are nonvoid. Moreover for any 
i G I and j G I', bij[x) = bji{x) = 0. Indeed, {DB{x)D~^)ij = didj^bij{x) and 
{D'B{x){D')~^)ij = d[{d'j)~^hij{x), yielding bij{x) = 0 and, similarly, bji{x) = 0. 
Hence, the system is decomposable. □ 


The proof of Theorem |3.5| shows that there is considerable choice for the deformed 
monodromy M^o, indeed an open set’s worth. We note that it provides a local 
normal form in the sense that given the choice of Stokes matrices , Cj", one then 
has an open set of suitable Mao to choose; once one has done this, the resulting 
Mao then determine the normal form (3.1) uniquely. One can use the 






diagonal automorphisms to further normalise the connection, by normalizing some 
off diagonal coefficients to one. 


3.2. The special case of an irreducible irregular singularity. The irreducibil- 
ity we have is for the system as a whole on CP^; on the other hand one has gener- 
ically that the irregular singularity itself is irreducible. Bolibruch showed the fol¬ 
lowing theorem 


Theorem 3.11. or Theorem 20. f of m) A linear system of differential equa¬ 
tions with a nonresonant locally irreducible singularity of Poincare rank k is lo¬ 
cally holomorphically equivalent to a polynomial system x^'^^y' = A(x) ■ y, where 

Mx) = ELo 


Corollary 3.12. A linear system of differential equations with a nonresonant lo¬ 
cally irreducible singularity of Poincare rank k is locally holomorphically equiva¬ 
lent to a unique normalized polynomial system x^^ y' = A{x) ■ y, where A{x) = 
with a normalization as in Lemma 


3.10 


Proof. We can of course suppose that Aq is diagonal. Moreover, an irreducible 
bundle with connection is in particular indecomposable. Hence, it is possible to 


apply a further normalization of (n— 1) off-diagonal coefficients using Lemma 3.10 


□ 


The condition of irreducibility in Bolibruch’s theorem in some sense explains 
the need for introducing the extra singularity in our normal form; we need some 
freedom so that the monodromy around oo be sufficiently transverse to the diagonal 
structure at the origin coming from >1(0). 

3.3. Deformiug iu a fixed DS domaiu; coutiuuity at e = 0 iu a disk Djj. 

We would now like to vary e in a fixed DS domain Sg. Since the DS domain is fixed. 
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we drop the index s. We first work on a disk in x-space. Later, in Section |3.4| 
we will extend to the whole of CP^. 

We begin by clarifying our notion of a continuous family of bundles plus connec¬ 
tion, as e varies. Indeed, if one has a family of bundles given in pictures A or B 
of subsection 2.8 above, one can imagine that there could be difficulties in under¬ 
standing continuity as one varies e, since the transition matrices are defined over 
sets varying with e. Within the DS domain Ss, this poses no problem, as the open 
sets are varying quite smoothly. But one should be a bit careful as one moves in to 
e = 0, or more generally, as one is going to the boundary of the DS domain, at a 
point where the singularities merge together, for example with two poles becoming 
a double pole. 

On the other hand, if one thinks of our pair of (bundle, flat connection) as 
a family of flat continuous connections in a global (not necessarily holomorphic) 
trivialization (let us call this, picture D), then it is fairly easy to see what a good 
notion of continuity should be. We first must fix the singularities of the connections; 
this amounts to choosing conjugacy classes of the singular part of the connection, 
up to holomorphic gauge transformation. This polar part must vary continuously 


For our example, in the disk Djj around the origin, we first choose a formal 
normal form, which must itself be deforming continuously in e: 


(3.3) 


An.e — 


Ao(e) -f ....Afc(e)x* 
Pe{x) 


dx 


Our continuous family Ac will have poles that are holomorphic conjugates of 
our normal form, and will have in addition a finite piece, living in an appropriate 
function space, with both (1,0) and (0,1) components: 


(3.4) 


Ac - f{x, e)An,efix, e) ^ = a^’° -k 


, 0.1 


with f{x,e) S Mat(n x n,C) invertible, holomorphic in x and varying continously 
in e, and with -|-a°’^ lying in a Sobolev space and varying continuously in 
e, as elements of the Sobolev space; note that this does not imply that the and 
aj-g themselves continuous. Thus, while the (1,0) (holomorphic) term of Ac 
has singularities, the (0,1) term is bounded and continuous, again as a family 
of elements a°’^(e) of with 1 < q < 2. 


Remark 3.13. Some comments are in order as to the choice of space W^ '^, a 
space of functions in L'^ with one L'^ derivative. One wants to solve for a gauge 
transformation, gf^dg^ = —which will take us to a holomorphic gauge, and 
we would like to have the solution be continuous, so that, for example applying 
the gauge transformation does not change the topology of the bundle; this suggests 
setting up our problem as finding an inverse image under a map —)■ 

since one has the Sobolev embedding theorem C C®, for q > 1; of course, over 
the subset where is smooth, elliptic regularity will tell us that g^ is also smooth. 


Remark 3.14. While the notion of continuity allows for the gauge freedom given by 
f in (3.41, in our case, however, this freedom can be normalized away by preparing 
Ac so that it has the same polar part as An,c, and we only consider the equation 
with f = 1: 


( 3 . 5 ) 


Ac — An,c = + a°’^. 
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Now consider a family of connections in picture A or B, i.e., given by Stokes 
data (A), or Stokes data shifted to the normal diagonal form (B). In picture B, 
explicitly, one first builds a bundle by glueing trivialised bundles on the sectors 
n Dij by symmetries of the normal form coming from the Stokes matrices on 
intersections; the connections are constructed by putting the normal form An,e on 
the sectors n D/j, and patching together with a partition of unity, obtaining a 
C°° connection on Djj; initially this is on the complement of the singularities, but 
the asymptotics of the glueing functions are such that the bundles extend to the 
singularities. More explicitly: 


Definition 3.15. (Definition of the pair (bundle, connexion Ae)) Let us describe 
the bundle, with the transition functions coming from the Stokes matrices; this is 
essentially picture B above. Let be the standard fundamental matrix solutions 


for the connection given by (3.31 on the which are analytic continuations 

of each other, except on Then, as in (2.131, (2.141, (2.15), the transition 

functions are of the form 


Ct = F+= dl 


(3.6) 




■N-, 






= +iV^, xGO^^„ 




=nfiUu 




Here, the transition functions of ( |3.6[ ) are all of the form Dg + iVg(x), with Dg 
diagonal, constant in x and analytic in e, and TVg either strictly upper or lower trian¬ 
gular and tending to zero at the singularity. Since + fVg {x) is an automorphism 
of Mn,£ near the origin, it satisfies 


(Dg + lVg(x))M„,g(Dg + N,{x))-^ + dN,{x){D, + iVg(a:))-i = A 


(where d is the derivative in the x direction only), whence 


dN^{x) - [An,e,Ne{x)] = 0, 

since [D^,An,e\ = 0. This equation involves the formal. Abelian connection, and 
decouples each entry of the matrix iV; in consequence, the entries i < i' ol 
N^{x) (to fix ideas, in the upper triangular case, with Re(t) —)■ —oo), an upper 
triangular matrix with zero diagonal terms, satisfy 


(3.7) 


ni^i>{x,e) ~ Cjy/ exp((Ai(a;£) - Xi>{xe))t ), 


where 

Ao(e) H-h Ak(e)x’^ = diag(Ai(x, e),..., A„(x, e)). 

The asymptotic decay would be an application of the theorem of Levinson given 
above (Theorem 2.5), but in the simpler scalar case, since the equations for each 
are decoupled. 

The important ingredient is that Re(Aj(xf) — \ji{xg)) > 0 for j < j'. Thus, 
Aj(xf) — \j'{x() = with G (~f i f )• Let us suppose that we consider 

t along a line of slope 6 G (~f ? f )> he t = with t' G ]R“ (since x^ is of a-type 

in the upper triangular case). Then 


exp((Aj(x^) - Xji{xe))t) = exp 


and 

(3.8) 


\nj^y{x,e)\ < C'(e)exp(-/3|t|), 
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for some /3 > 0 provided that |0| < <5 with 5 > 0 sufficiently small so that + 
9 € (—f 7 f )• It suffices to take /3 = {rjji cos{(j)jj> + 9)). (The number S 

appearing here was the limit of the slope for the strips in f-space in |^.) A similar 
estimate holds for the entries of N^{x) in the lower diagonal case when Re(t) > 0. 
This gives transition functions which tend to constant diagonal matrices in the 
limit, and so the bundle extends to the punctures. 

We now can build the connection. On each sector, we have the normal form 
connection An,e- Let /ig be a smooth function, which is zero on \ 1 on 

\ Oq,, and takes values in the interval [0,1] on the intersection. We then set the 
connection on 0^ U 0,3 to be given by the gauge transformation on the overlaps: 

( (£>£ + heix)N^{x))An,eiD<^ + heix)N^{x))~'^ 

(3.9) Ae = < +d{he{x)N^{x)){D^ + he{x)N^{x))~^, on O^ 

£, on 0,3 \ Oq,. 

Note that A^ = An,e on (0(j\0^)U (0^\0 q,), and that the connection is everywhere 
flat. Iterating the process for all intersection sectors allows defining Ae globally on 
D^j, with the caveat that there might be some singularities in the connection at the 
singular points, a problem we will discuss below. 


Proposition 3.16. Let the Stokes data vary analytically in the intersec¬ 

tion of the DS domain Ss with the polydisk Dp in e-space, with a continuous limit at 
the boundary of the DS domain, near the points at which some of the singularities 
coincide (i.e. A(e) = 0/ Then, choosing the functions appropriately, the pair 
of Definition 3.15 defines a continuous family of pairs (bundle, connection) on D/j 
in the sense given above, that is as a map of Ss H Dp into 

Proof. For the visualization of the sectors in the t-coordinate we refer to Figures]^ 

El HKD and [8| 

Our Stokes data is defined on a family of bundles over the disk D^ in C. One now 
wants to check that the result, when converted to “picture Z?”, gives a continuous 
family. This is done by passing from trivializations defined on our various open 
sets to a common trivialization, using smooth but non analytic functions on the 
overlaps of the open sets. As noted above, within the DS domain, everything (the 
Stokes matrices, and the open sets themselves) is varying smoothly, and obtaining 
continuity poses no problem, and indeed is quite classical. 

The problems occur when one approaches the discriminantal locus A(e) = 0, 
and indeed, only then in a neighbourhood of the singular points. Note that A = 0 


is invariant under (2.6l. Near a regular point of {A = 0}, we can reparameterize 
e I— >■ (77, 77 '), with 77 = A. 

We have glued all our local definitions together in a single patch as described 
in Definition 3.15 One then considers Ae — An.e, referring to (|3.9|. We first want 
to estimate its L'^-norm. Since {De + he{x)Ne{x))~^ is of the order of 1, this is 
tantamount to estimating the norm of {Ae — An,e){De + he{x)Ne{x)), which is 
equal to 


{De + he{x)Ne{x))An,e “ AnA^e + he{x)Ne{x)) - d{he{x)Ne{x)) 
= [he{x)Ne{x),An,e] - d{he{x)Ne{x)), 
and so, given that dNe{x) + [An,e, ^"£( 2 ;)] = 0 , the norm of the quantity 

Ne{x)d{he{x)). 
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Thus the behaviour of our chosen he as we approach the boundary is quite crucial. 
This is where the t-uniformisation of the plane comes into play. Under this uni- 
formisation, the sectors have a horizontal part which, together with the width 
of the sector, becomes large as e tends to zero. We can of course manage that the 
intersection of two adjacent sectors be a strip V of uniform vertical width 1. When 
we approach a boundary point e' of Ss belonging to {A = 0}, the same occurs: 
some sectors (not all, only the ones attached to a multiple point) become large 
as well as the horizontal part when e —>■ e'. Hence, we can use a function in the t 
plane, which just depends on the imaginary part of t, and goes smoothly from zero 
to one from one side of the intersection strip V to the other, so that dh^ is supported 
on V. We then have that dhe is of the order of one, and we want to estimate the 
norm Ne{x) on the strip V of vertical width one, going to infinity. This is, however, 
already done; from dNe{x) + \An,e, fVe(a:)] = 0 , as we noted above, but now passing 
to the f-parametrization, the entries of W are of order exp(—Q;|Re(t)|), as noted in 
(3.8l. One then finds, taking into account the change of variable, that the quantity 


one wants to bound for the L'^-norm is then 


(3.10) 


C(e) / exp{—qa\t\)\pe{x)\'^ ^dtAdt, 


which indeed remains bounded, independently of e, if 1 < g < 2 , as |pe(a::)| is 
bounded near t = oo. Similarly, for the L'^-norm of the derivative, one has that 
-^{dhe) is also bounded, and so the derivative ^ is of the order of p~^. The 
quantity one wants is then 


(3.11) 


C{e) / exp{—qa\t\)\pe{x)\^ '^'^dtAdt, 


which requires a bit more work to bound. To fix ideas, we can consider the lin¬ 
earized equation for x(t) at a singular (say attractive) point, which is i; = p'^{xi)x, 
giving a solution x{t) = exp(p'(a:^)^), with Yie{p'^(xi)) < 0 , and so in the linear 
approximation |pe(a:)| ~ \p'^(xi)\exp(p'^(xi)t)^ and the integrand one is looking to 
bound is \p'^{xt)\^~‘^^ ■ exp((—qa -|- {2q — 2)p'^{xt))\t\). Now, as e gets small, so do 
the derivatives p'^ixi) at the zeroes of p, and so for e small we get an exponentially 
decreasing integrand. Likewise, at a multiple zero, the corresponding approxima¬ 
tion gives x{t) = A, for some 7 and so the integral also converges. This argument, 
of course, is only heuristic; for example, there is no bound on the constant C'(e) as e 
approaches the discriminant locus. Nevertheless, it indicates that some estimate is 
possible, and indeed, it is provided by the following lemma, which gives a uniform 
bound in e over open sets. 


Lemma 3.17. Let X£{e) be a zero ofpe{x). Suppose that there exist a, c > 0 such 
that 

Re ( and |pe(a:)| > c|a; - a;^(e)|'‘+^ 

for \x — xi{e)\ < S, and e € D. Let x{t, e) satisfy 

=Peix{t,e)), 

and let it converge to xg(e) as Re(t) -A -too in a strip. Then there exist constants 
K, K' > 0, with \x{t, e) — xe{€)\ > K exp(—at) and Pe{x{t, e)) > K' exp(—(fc -I- l)at) 
for t > t(e). 

Similar estimates apply in the case o/Re -A — 00 . 
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Proof. The first constraint gives 

- .,))) = ite ((7^^) > 

and so, integrating, Re(ln(x — X{)) > C — at, and jcc — xi\ > K exp(—at), whence 
\pe{x{t,e))\ > K'exp{—lk + l)at). □ 


Returning to the theorem, for e sufficiently small one can make the bound a 
of the lemma sufficiently small for the estimate (3.111 to converge. We then note 
that, on finite strips, all of our constructions are continuous in e; this, plus the 
exponential convergence, gives us a continuous family in e in the function space 

□ 


From this continuous family, in (picture D) on Ss, one can pass to a 

holomorphic description, in the neighbourhood of the origin, of the connection in 
a single holomorphic trivialization. One can proceed as in the paper of Atiyah and 
Bott jT], p.555. We note that starting from picture D, the aim is to find a gauge 
transformation solving 

(3.12) gf^dg^. = -a°’^ 


We know that there exists a gauge transformation u to a holomorphic gauge for 
e = 0. Gauge transforming the whole family with u, we can assume that we are 
deforming from a holomorphic trivialization at e = 0 , so that a®’^|e=o = 0 ; in short, 
we are solving, as e varies in a small set, for in a neighbourhood of the origin 
in Viewing 

(3.13) 5 , gf^dg^ 


as a map of Sobolev spaces —)■ on a disk around the origin, one wants to 

appeal to the implicit or inverse function theorem on Banach spaces. To do this, 
we compactify, so that the map on function spaces over CP^ is Fredholm. We first 
extend the bundle trivially over CP^. We can of course suppose that e is sufficiently 
small so that the singularities all lie within {|a;| < ^}. We use a C°° function with 
bounded derivative 


f{\x\) 


1 , |a;| < f, 
0, |a:| > R. 


and we extend Ae and An,e to f{\x\)Ae and f{\x\)An.e respectively, on a family of 
trivial bundles {E^} on CP^. 

The linearization at 5 = Id is the family in e of elliptic Dolbeault complexes 

I—)■ dg^, mapping the sections of a trivial bundle to the sections of the tensor 
product of the same trivial bundle with the line bundle of (0,1) forms. Its kernel 
is the family of global holomorphic sections 7F°(CP^,End(Ae)) = gl{n,C), and 
its cokernel is the first Dolbeault cohomology group iJ^(CP^,End(ife)) = 0 (see, 
e.g. [ 8 ], chap. 0), and so the Fredholm map (3.131 is locally surjective near the 
identity as a map from the Sobolev space of sections of Aut(ifg) with two L‘^ 
derivatives to the space W^’’^ of sections of (0, l)-forms with values in End(£'e) and 
one derivative (here <7 > 1). Hence, the map (3.131 is surjective, and its kernel 
is given by the constant sections. Asking for example that gig be orthogonal to the 
constants gives, by the inverse function theorem on Banach spaces, (e.g. Lang |15| . 
p. 15) a unique solution gg for (3.121 restricted to a suitably small open set in 
e-space. The details follow Atiyah and Bott [T]. 
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Transforming Ae with we obtain a family of connections Ve = ^ + 

dg^ g~^. The (0,1) part of Ve, namely {g^a^’^ + dge)g^^, vanishes by the construc¬ 


tion of g^ as solution of (3.12). Since Ve has been obtained from a flat connection 


by gauge transformations, it is flat. Its flatness then ensures that it is holomorphic. 

Hence, we obtain a family of connections on a disk depending analytically on 
e G Ss with continuous limit at points of the closure of lying on {A = 0} (this 
includes e = 0). 

3.4. Connections on CP^ over a DS domain. The previous section gave us a 
family of bundles in picture C over the product of a disk D/j in C and a closure 
of Ss n Dp, with a singular connection in the disk direction. (Though we used an 
extension to CP^ in the previous lemma, this is not the extension we will use here. 
) 

Indeed, we proceed as for the bundle at e = 0 and glue in two Fuchsian singulari¬ 
ties; the result will be a deformation of our normal form for e = 0. The deformation 
will then also be irreducible, and the underlying bundle will be trivial. The normal¬ 
isation extends to these deformations also. Writing out the connection in a global 
trivialization, we will then define the analytic normal form for our singular rank k 
systems; we suppose th at w e have already produced our extension of the connection 
to CP^ at e = 0, as in (3.5), and so have a monodromy at infinity and residue 
at infinity A^: 


Theorem 3.18. Suppose given 

• a DS domain Ss, 

• formal invariants given by diagonal matrices Ao(e),... Afc(e) depending an¬ 
alytically on e in the intersection of a polydisk Dp with the closure of Ss, 
such that Ao(0) has distinct eigenvalues satisfying o, determining a for¬ 
mal normal form (1.8), 


collections of normalized invertible (Stokes) upper (resp. lower) triangular 
matrices (resp. j = depending analytically on e G 


(3.14) 


Dp n S's, with continuous limit at the boundary points, and determining a 
monodromy 

M{e) = . 

a generic matrix M^o in a neighbourhood of the identity with distinct eigen¬ 
values, representing a conjugacy class with distinct eigenvalues, with all en¬ 
tries nonzero, close to the identity, and such that Mqo has distinct 

eigenvalues, 
conjugates 

Moo(e) =r(e)Moor(e)-i 

with non-zero entries for e in a small neighbourhood of 0 for some suitable 
analytic function T : 5's — >■ G(n, C) with continuous limit at the boundary 
points, and such that r(0) = id. 

Then, restricting p if necessary, there exists a family of irreducible rational linear 
differential systems 

Ao(e) + Ai(e)a;Afe(e)a;^ A[e) 


(3.15) 


(3.16) 


y = 


Pe{x) 


X — R 


■ y = B{x,e) ■ y, 
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depending analytically on e G 
(3.1 1 for e = 0) , and with 


Ss with continuous limit at A(e) — 0 (of the form 


formal normal form at the origin (1.8 1 , defined for e G S's H 


generalized Stokes matrices which are given by the matrices and 
j = l,...,fc, 

and monodromies Mule), M^{e) around such that Mji{e)M^(e) = 


M{e), with M{e) defined in (3.141. 


The automorphisms of the system are multiples of the identity. Furthermore, 
acting by an automorphism of the formal normal form (a diagonal matrix, constant 
in x), which changes both the Stokes data and our rational system, it is possible 
to normalize the coefficients of n — 1 non diagonal entries of suitable off-diagonal 


entries of B(x,e) to 1, thus leading to a unique normalization of (3.1 1 for each 
equivalence class of the Stokes data under automorphisms of the formal normal 
form. 

Proof. We have already constructed a family of bundles with connection (ife, Ve) 
over Dij, which is a continuous deformation of (i?o,Vo); we glue to this, as for 
e = 0, the bundle on the disk Doo = {|a;| > R/2} with two Fuchsian singularities at 
R, oo with monodromies. This is done by choosing a continuous family of bundles 
over Doo x (Ss n Dp), trivialized at a: = 3i?/4, with a connection with Fuchsian 
singularities at x = i?, oo, monodromy M{e) around the circle of radius 3i?/4 
defined in (3.14), and monodromy around infinity given by Moo(e). The bundles 
with connection over D^j and Dqo are then glued in the standard way, starting from 
the base point. This gives the desired global bundle with connection, deforming the 
case e = 0. Since a small deformation of a trivial bundle is trivial, we can pass to 


a global trivialization. We then get a connection which has the form (3.16) above. 


As noted, for this connection, after diagonalizing the leading term, we can then 
normalize to 1 the same (n — 1 ) non diagonal terms as in the case e = 0 . □ 


Our Stokes matrices for a connection with Fuchsian singularities depend on the 
DS domain, as the matrices depend on the way that the singular points are tied to 
infinity via the separatrices; however, there is one invariant that does not depend 
on this data, and that is the monodromy representation. The importance of the 
monodromy representation is that it essentially determines a connection with poles 
of order one, up to some discrete choices of the polar parts of the connection: 
basically, the eigenvalues of the monodromy determine the eigenvalues of the residue 
matrices at the Fuchsian singularities up to integers. 

Proposition 3.19. Suppose that two pairs (trivialized bundle, connection with 
Fuchsian singularities) on CP^ have, in our context 

• the same singularities, given by the zeroes of Pf^{x),as well as R, oo; 

• the same formal invariants Ao(e),... Afe(e); 

• the same conjugacy class of residues at R and oo, given by the classes of 
Aii{e), Aao{e), with distinct eigenvalues; 

• fixing the base point, the same monodromy representations, (the same, not 
simply conjugate), with distinct eigenvalues around R, oo. (Note here that 
the trivializations on the fiber above the base point are given by the trivialized 
bundles, and so the monodromy map is uniquely defined.) 

Then they are isomorphic. 
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We emphasize that our bundles are trivialized, not simply trivial. The proof is 
given in the course of the proof of Theorem |4.6[ 

4. The compatibility condition: glueing DS domains. 

Given our Stokes matrices g , g , with same limit at e = 0, independently 
of s, we want to realize a corresponding differential equation (with a fixed formal 
normal form) over a polydisk in parameter space. We have already done this over 
our DS domains, and want to sew the results together. This involves three steps: 

(1) Glueing over DS domains. If the Stokes matrices are arbitrary on each DS 
domain, then there is no reason why the realized families over DS domains 
Ss and Ss' should be analytically equivalent one to the other over Ss H Sg' ■ 
A suitable compatibility condition is necessary to ensure that this is the 
case and so to allow a glueing of the different DS realizations in a uniform 
family depending on e € Eqi where Eq is the complement of {A = 0}; 

(2) Showing that this extends to the generic locus of A = 0, where A is the 
discriminant of Pe(a:); this generic locus is the set of e for which Pe{x) has 
one double zero and the remaining roots are simple; 

(3) Extending to the rest of the polydisk by appealing to Hartogs’ theorem. 
This section is concerned with the first step. 


Recall that Eq is covered by the Ck DS domains: we define the compatibility 
condition on the intersection of the DS domains. So far, given our Stokes matrices 
^pe,Ss > e S, > realized our system abstractly as a pair (bundle, 

connection) on the Riemann sphere, and showed that we could then realize it also 
as a singular connection on a globally trivialized bundle. For this connection, we 
had some freedom on the choice of the monodromy M^(e) around infinity defined 
in (3.15) and the monodromy around R is determined accordingly. 

No matter how it is presented, such a bundle over CP^, equipped with a singular 
connection on the complement of the zeroes xi(e), ..,a:fc+i(e) of Pe, and of Xk +2 = 
R,Xk +3 = oo, comes with a monodromy representation 

Afe : 7ri(CP^ \ {a:i(e), ..,Xk+ 3 }) GL(n,C), 

defined up to global conjugation: one chooses a base point, then integrates. We 
note that since Moo(e) is given in (3.15) and M(e) is given in (3.14), one is in 
essence looking at the monodromy of the restriction of the connection to a disk 
Dsh of radius ^ : 

Me ■ 7ri(D3H \ {a:i(e), ..,Xk+i{e)}) ^ GL{n,C). 

4 

We do note, however, that the automorphisms of the bundle on Dsh a priori 

4 

just form a subgroup of the diagonal matrices and not necessarily scalars; it is only 
when one goes to CP^ that one is forced to an indecomposable representation. 

Also note that as one changes DS domains, the open sets on which the transition 
matrices C^^e,Ss ’ e Ss defined change; a bifurcation occurs. If one 

computes the monodromy along a loop, one gets an ordered product of inverses of 
transition matrices of the open sets that the loop intersects, taken in the inverse 
order of intersection (see for instance Example 4.7 below). In a different DS do¬ 
main, for the same loop, one gets a different product, even though the monodromy 
representations must be the same. 
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Figure 9. A self-intersection of Sg providing a tubular neighbor¬ 
hood of A = 0. 


Compatibility Condition 4.1. Let us assume given our data of a bundle plus 
connection, in picture A (or B), over our family of DS domains Sg- The Compat- 
iblity Condition is as follows: 

• We fix a base point x = ^R, and choose identifications of the fibers of our 

bundle over that base point, holomorphically and continuously up to the 
boundary in e G SgH Sgi. Then, for each e in SgH Sgi, we ask that the mon- 
odromy representations Ale, At' defined by Cl^^e s,) ^/"e s, j e s 

s , j s / ’ e s / equivalent, that is conjugate by an invertible 

matrix Ge = Ge(s, s') depending only on e, holomorphically, continuous up 
to the boundary points of each DS domain, and such that Go(s, s') = id. 

(4.1) TWe.s, = GgM ,, s ,, G :\ eGSgHSg ,. 

• The Ge(s, s') form a cocycle (Ge{s, s')Gg{s', s") = Ge(s,s")) and this cocy¬ 
cle is trivial: there exist invertible matrices re(s) depending analytically on 
e € Sg with continuous limit at the boundary points of the DS domain such 
that 

(4.2) Ge(s,s')=re(s)-lre(s'), 
and ro(s) = id. 

• The Sg forming a covering on the etale sense, condition also holds when 
Sg is a ramified DS domain as in Figure and we consider a connected 
component of its self-intersection (this happens for instance when k = 1, 
and also when considering neighbourhoods of regular points of {A = 0}/ 
In that case, one also asks that the matrices be a trivial cocycle in the 
sense that the matrix Ge(s,s) be equal to the product re(s)re(s)“^, where 
r,r represent the two different branches of the function re(s) as one turns 
around the divisor. 


Remark 4.2. The Compatibility Condition \4.1\ is necessary. Indeed, when con¬ 


sidering an analytic family (1.31, it is obvious that the monodromy representa¬ 
tions conjugate since they are two representations of the 

monodromy of the same family (1.3l. For the same reason, the cocycle of the 
Gg{s,s') is trivial. What needs to be proved is the limit properties o/Ge(s,s') 
and re(s). Let us suppose for instance that the given set of loops starts at the 
base point xt = ^ in ^ ^ n flies'- shown in [3] that well chosen nor¬ 

malizing changes of coordinates H^e s (fibered in x!) over ^fgg transforming the 
normal form (1.81 into the family in have a limit at points of {A = 0}. We 
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can of course consider that the monodromy representation Me,Ss expres¬ 

sion of the monodromy maps in the canonical basis over the base point X},. This 
canonical basis is transformed into a basis given by the columns of for 

(1.31 and A^e,s, is the monodromy representation in that basis. We want to write 
the monodromy Me,Ss ® fixed basis, independent of e and Sg- We choose as 
a fixed basis the one given by the columns of Hff^ixb). Then the monodromy 
is of the form rg(s)A^j 5 ^ (r£(s))“^, where the columns of the matrix (rg(s))“^ 
are the column vectors of Hj^Q^Xb) written in basis formed by the column vec¬ 
tors of ^ g (xb), namely H^Q{xb) = Hjf ^ ^{xb){Tg{s))~^, which yields obviously 
r£(s) = [Hff Q{xb))~^^ s(^b), from which the result follows. 


Remark 4.3. Since we are deforming from the same connection at e = 0, one 
could have dropped the condition that the matrices Gg{s, s'),Tg{s) take value 1 at 
e = 0. We would then have allowed the values at e = 0 to lie in the automorphisms 
of the formal normal form. Then one does not necessarily have the same Stokes 
data, but equivalent Stokes data at e = 0. 


Lemma 4.4. The Compatibility Condition \4-l\ allows the monodromy representa¬ 
tions to be made equal (instead of just equal up to conjugacy), in a consistent way. 
Indeed, acting by Fg in the fibers over our base point, and so setting 

Mg,s.=Tg{s)Mg,sr:\s), 

one obtains that the transformed G'g(s,s') are the identity. 


Of course, once one has done this, the monodromy on a loop is no longer an 
ordered product of Stokes matrices and gate matrices and their inverses, but rather 
a conjugate by Fg of this product. 

We can extend the normalised representations M.e,Ss to the full complement of 
fc + 3 points in the Riemann sphere; one has fairly immediately: 


Proposition 4.5. Let Sg and Sgi be two intersecting DS domains. The Compatibil¬ 


ity Condition j.l implies that the full monodromy representations Afe,Ss ^^id Afe^s 1 
around the fc + 3 singular points are the same (or their conjugates Afe,Ss bind Afe^s / 
by Fg(s) and Fg(s') respectively), provided we choose the monodromy at 00 as in 
(3.151 (with Mao having all nonzero entries) in the realization over each DS domain 


given by Theorem \3.1^ 

Proof. In other words, if we conjugate by the Fg(s),Fg(s') (change trivializations 
at the base point), so that the monodromy representations are the same, we then 
add in the same monodromy matrix at infinity. More explicitly, we consider the 
monodromy Afs^(e) and AAs^iie). By the Compatibility Condition 
that 


4.1 


we have 


Mg.s, =Fg(s)Ms,(e)(Fg(s))-i =Fg(s')Xs,,(e)(rg(s'))-^ 

When constructing the realizations, the choice of monodromy we make at infinity 


in (3.151 with the function Fg(s) defined in (4.2) guarantees that the monodromies 
around infinity in the glued domain {|a::| > U { 00 } calculated from the same 
sector LIT are the same, indeed constant and equal to M^o- Then equality also 


follows for the monodromies Mfi and M^ s , around x = R. 


□ 
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Theorem 4.6. Let 8$^ and £$ , be two normalized realizations over CP^ as con¬ 


structed in Theorem 3.18 If the Compatibility Condition 4-1 is satisfied, then 


and £s , are equal on Sg H Sg 


Proof. We have seen that we could consider that the monodromy representations 
of and £s , are equal. We consider a base point X}, near infinity and loops 
7 ^ G ni(CP^,a;f,) surrounding xg, alone in the positive direction, i = 1,..., fc + 3. 
We can suppose that the points are numbered so that 7 f ^ ■ - 1^+2 homotopic 
to 7 fc+ 3 . Let Xg (resp. Xgi) be a fundamental matrix solution whose columns are 
eigensolutions at oo of £s„ (resp. £s„,)- 

Let us call the monodromy of Xg and Xgi along 7 ^. 

One first has that the two realisations are isomorphic on the complement of the 
singular set. This basically is simply the fact that the induced monodromy on 
the bundle Hom(5g^,£s /) 3 .cts trivially on the analytic continuation of the identity 
map. More explicitly, we consider the map y 1 —> P{x)y = XgiX^^y. We need to 
show that P{x) is well defined. The analytic extension of Xg (resp along 7 ^ is 
given by XgMg^xt (resp. XgiMg^fi). Then P{x) is well defined since 

{Xg,Mg,,,,) {MfilXf^)=Xg,Xf\ 

The next step is to show that P can be extended analytically to the singular 
points. Let us start by showing that P can be extended at Xk +3 = 00 . Our 
singularity is non resonant, i.e. no two eigenvalues of the residue matrix differ 
by an integer, thus ensuring that the monodromy has no multiple eigenvalues and 
is hence diagonalizable; it is a generic situation for Fuchsian systems, in which 
the monodromy around the singularity is determined by the eigenvalues of the 
residue matrix and the system decouples as a sum of one-dimensional systems 
( 0 ). Thus, if sxe the residue eigenvalues at Xoo, then the columns 

of Xg (resp. Xgi) which are eigensolutions are of the form Wgj{x) = x~^^ggj{x) 
(resp. Wg'j{x) = x~^iggij{x)), where ggj (resp. gg'j) is analytic and nonzero 
in a neighborhood of Xoo- Moreover, the matrix Lg (resp. Lgi) with columns 
gsp{x),... ,gs^n{x) (resp. gg>p{x ),..., (/s',n(a:)) Las nonzero determinant for x close 
to Xoo- We have that P{x)wsj = Wg/j. 

Then 

iXg = Lsdiag(a;“^S...,a;“^"), 
yXgi = Lgi diag(a;“'^i,.. 

Hence, P{x) = XgiX~^ = LgiLj^ is a nice analytic matrix with nonzero determi¬ 
nant in the neighborhood oi x = Xao- 

Let us now consider a singular point xg at which the monodromy Mx^ is diag¬ 
onalizable (this is in particular the case for xg — R from our construction), and 
let N G GL{n,C) be such that N~^MxiN = D, where D is diagonal; we place 
ourselves again in the generic situation at which there is no resonance (i.e. no 
two residue eigenvalues differ by an integer). The matrices XgN and XgiN are 
still fundamental matrix solutions of £g and £gi, and their columns are eigensolu¬ 
tions for the monodromy around xg. If /ri,...,/i„ are the residue eigenvalues at 
Xg, then the eigensolutions are of the form Wgj = (x — xg)^^gsj{x) (resp. Wg'j = 
(cc — xg)^^ggij{x) ), where ggj (resp. gg'j) is analytic and nonzero in a neighbor¬ 
hood of Xg. Moreover, the matrix Qg (resp. Qgi) with columns gsp{x),... ,gg^n{x) 
(resp. ggi^i{x),... ,gg\n{x)) has nonzero determinant for x close to xg. This comes 
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from the fact that each system is diagonalizable near Xi and the eigensolutions for 
the diagonal form are given by C{x — X()^^ej, where C G C* and Cj is the j-th 
vector of the standard basis. Then 

fXgiV = Qs diag((x - xiY '^,..., (a; - 
yXs'N = Qs' diag((a: - xeY\ ..., (a; - x^YY- 

Hence, P{x) = Xg'Xj^ = Qs'Qj^ is a nice analytic matrix with nonzero determi¬ 
nant in the neighborhood of a: = 

We have built an equivalence depending analytically on (x, e) (since it is the case 
for Xs and Xg') on the domain CP^ x S'\ {(x, e) : Pe(x) = 0, e resonant}. From our 
hypothesis, the set of resonant values of e is of codimension 1 in e-space. Hence, 
the set {(a;,e) : Pe(x) = 0 ,e resonant} is of codimension 2, and we can extend the 
equivalence to it by Hartogs’ theorem. □ 

4.1. Generic bifurcations and the compatibility condition. Now suppose 
that e lies in a connected component Sg H Sg' of two DS domains: the bifurcation 
from Sg to Sg! is precisely obtained by changing some angle 9 with which one goes 
out to infinity; this can switch some singular point (s) Xf, from a-type to tu-type 
(or the converse), or can change the points of attachment of the separatrices at 
infinity to the singular points. We would like, in this subsection, to illustrate how 
this impacts on the generalised Stokes matrices. 

One thus has a bifurcation; in codimension one (the generic bifurcation) there 
are basically two types of bifurcations in e that force a change in DS domain. Both 
involve going through a homoclinic connection between two of the separatrices at 
infinity. In the process, one gate sector is untied from its endpoints and each end 
is then tied to a new endpoint. 

• Outer connection: a zero of Pe(x) changes from being a-type to w-type; 
this occurs at one end of a branch of the skeleton as in Figure 

• Inner connection: All the zeroes keep the same a-type or w-type as before, 
but the attachment of the separatrices emerging from infinity to the zeroes 
of Pe changes as in Figure pT] 


Example 4.7. An example of the Compatibility Condition. The Compat¬ 
ibility Condition means the following: given any set of loops 'yxi starting from a 
base point and going around xi in the positive direction as in Figure \1^ we can 
write the abstract monodromy around that loop as a product of matrices 

g , C^^ g , e s their inverses corresponding to the intersection sectors 

that are crossed around that loop. The set of matrices monodromy 

representation compare with / • This gives for the example of 


Figure 12 


= (C,\) 


-irG 




fCl 


M. 


X2,Ss 


M. 


X2,Ss 


M. 


X3,Ss 


M. 


^3,Ss 


_ (r'G 'i —\—^nG r^u 



Moduli space for generic unfolded linear differential systems 


39 





Figure 10. The two ways of passing through an outer homoclinic connection. 


4.2. Extending the realization to the generic locus of A = 0. The generic 
locus of A = 0 correspondes to two zeroes of coming together. For arbitrary k 
this reproduces a parametrized version of the case k = 1 studied in na. 

As one goes to this generic locus of A = 0, we note that the fundamental group 
of the complement of the zero locus of changes: it loses one generator. If one 
considers the representations defined by ^ 5 ^, there is indeed 

an issue, caused in particular by one gate element, say e s ’ gate is 

closing, or rather shrinking to zero. In general, also ^ g has no limit as one 

goes to A = 0. However, as shown in Sections |3.3| and |3.4[ the passage to “picture 
B”, with our choice of normalisation tames the C® / n o , to o = Id; if in 

addition the Stokes matrices g ^, g^ (or their modifications g ^, g ^) 
behave continuously when e tends to the divisor {A(e) = 0}, then we have a contin¬ 
uous family of connections, and our argument about a uniform passage to “picture 
C” goes through, giving a continuous limit. 
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Figure 11. The two ways of passing through an inner homoclinic connection. 


The same continuity is enforced when one has Ss self-intersecting (one is dealing 
with an etale covering) as one moves around the regular part of the divisor, as in 
Figure]^ for sectors as in Figure 13 We obtain a uniform differential equation on 
the tubular neighbourhood, with a limit at the core, as long as the compatibility 
condition for the monodromy representation is satisfied. 

We thus have a continuous limit at A = 0, which then must be a holomorphic 
limit at A = 0. 

On the rest of A = 0, we are in codimension two, and an appeal to Hartogs’ 
theorem suffices for the extension. 


5. The moduli space 

5.1. The realization theorem. 

Theorem 5.1. Suppose fixed 

• an integer k > 1, 

• formal invariants given hy diagonal matrices Ao(e),..., Afc(e) depending an¬ 
alytically on e in a poly disk Dp such that Ao(0) has distinct eigenvalues 


satisfying (l.ll. 


Now, suppose given for each DS domain Ss of radius p, collections of normalized 


invertible (Stokes) upper (resp. lower) triangular matrices 


('resp- 
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Figure 12. Two sets of sectors when k = 2 with a transition given 
by a homoclinic loop through the fourth quadrant and the corre¬ 
sponding monodromy groups. The permutation a is the identity 
(resp. the transposition on 1,2) for Ss (resp. Ss')- 



Figure 13. The sectors for a self-intersecting Sg in a tubular 
neighborhood of A = 0. 


j = 1,..., k, depending analytically on e € Sg with continuous limit at e = 0, inde¬ 
pendent of s, and continuous limit at generic points of {A = 0}. Moreover, suppose 
that the Compatibility Condition \4-. l\ is satisfied. Then, there exists an analytic fam¬ 
ily of rational linear differential systems (1.3 I for e € Dp, with formal normal form 











42 


J. Hurtubise and C. Rousseau 


(1.8 1 , and with Stokes matrices C^^acind j = over Sg- A par¬ 

ticular analytic family of the form (3.16) exists, with two extra Fuchsian singular 


points at R,oo, and the family can he uniquely normalised as in Lemma \3.10[ 

In the particular case where the system is irreducible for e = 0, a second realiza¬ 
tion exists of the simpler form 


(5.1) 


Pe(.x)y' = (Ao(e) + Bi{e)x H-h Bt,{e)x^) 


y, 


Proof. Let us call Eg C Dp the locus of in e-space of the polynomials with either 
distinct zeroes or one double zero and the rest distinct; it is the union of Eg (distinct 


zeroes; defined in (2.2)) and the generic points of A = 0 in parameter space (one 


double zero); its complement has codimension two. We have built an open covering 
{Va} of Eg, and on each open set Va, a realization by a linear system depending 
analytically on e € Va- The Compatibility Condition |4.1| ensures that for e € 
VanVp, the systems Ea and £p are analytically equivalent, and, once normalised, are 
the same. Hartogs’ theorem allows to extend the family of rational linear differential 
systems (1.3) to a uniform family over Dp; one uses implicitly the fact that a small 


deformation of a trivial bundle is trivial, and our normalisations Lemma |3.10| for 
uniqueness. 

The same type of arguments can be used verbatim in the particular case of 


an irreducible system at e = 0 and (5.1), since everything relies on a uniquely 
normalized realization at e = 0. Here we use Theorem 13.111 instead of our normal 
form. □ 


5.2. The moduli space. In the literature the “moduli space” is a universal space, 
typically finite dimensional, with all families of a given type of object given by 
mapping into that space, and often (for a fine moduli space) obtaining the family 
by pulling back. 

One could, of course describe the families directly, and in some sense that is 
what we are doing here. Our families of deformations will be, in effect, given by 
compatible families of Stokes matrices, and so by maps into the matrix groups. 
In that sense, we speak of the moduli space, even though what we are describ¬ 
ing, is, a priori, an infinite dimensional family of maps ^ g (e),..., ^ g (e) and 

C'j>o),e,s,(e)- These maps come with equivalences arising from actions of the group 
'Dn of invertible diagonal n x n matrices, corresponding to varying the trivializa- 
tions compatible with the flags on each sector, for each DS domain. Once one has 
normalized, this action reduces to one of the form 


(5.2) 


(-<L 



Kg{e)-^ 

/^U 



Kg{e)-^ 

riG 


f<G 



for Ks{e) a diagonal invertible matrix. We introduce the following equivalence 
relation on the collections of Stokes matrices: 


Definition 5.2. Two collections of (normalized) Stokes matrices {(7^^^ 

on a give DS domain Sg are equivalent if there exists a 
family of invertible diagonal matrices Kg(e), depending analytically on e G Sg with 
continuous invertible limit at e = 0 and at generic points of {A = 0} such that the 
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action (5.21 on the first collection of Stokes matrices gives the second. We denote 


the equivalence class by jcf ^ s, > • • ■ ^ ^ke s, > | 


eGS, 


Theorem 5.3. The moduli space under analytic equivalence for germs of generic 
unfoldings of nonresonant linear differential systems with an irregular singularity of 
finite nonzero Poincare rank at the origin and diagonal matrix Ao(0) with distinct 
eigenvalues satisfying ( 0 , is given by the set of tuples 


k, Ao(e),..., Afe(e), ■ ■ ■, 


eGS, 




s=l 


where 






k > 1 is an integer; 

Ao(e),..., Afc(e) are formal invariants given by germs of analytic diagonal 
matrices; 


for each DS domain Ss, 


{ 





are collections of equiv¬ 


alence classes of germs of invertible normalized (Stokes) upper (resp. lower) 
triangular matrices g , (resp. g ), j = 1,... ,k, depending analyti¬ 
cally on e G Ss with continuous limit at e = 0 independent of e, continuous 
limit at generic points of {A = 0}, and satisfying the Compatibility Condi- 
tion\4.1\ 


6. Ramifications 


Theorem 6.1. We consider a germ of family (1.3). If there exists a permutation 
matrix P such that the permuted Stokes matrices RC'jj.g P~^ have a common 
block diagonal structure with blocks of size ni,..., Um, n-i + • • • + Um = n for all 
j = 1,... ,k, for all Sg and for all f G {L, U}, then the germ of family is analytically 
equivalent to a direct product of germs of m families of linear differential equations 
on for each i. 

Proof. The modulus can be decomposed as a direct product of m moduli. We realize 
m families on C"b.. .C"™, having as Stokes matrices the corresponding blocks of 
size Hi and corresponding formal invariants. The direct product of these families 
has the same modulus as the original family. Hence, it is analytically equivalent to 
it. □ 


An immediate consequence of the realization theorem is the following normal 
form theorem. 


Theorem 6.2. A germ of family of linear differential systems unfolding an irregular 
non resonant singular point of Poincare rank k is analytically equivalent for an 
arbitrary choice of a non-zero R (independent of e) to a rational form 


( 6 . 1 ) 


Pe{x)y' 


Ao(e) + Bi{e)x + • • • + Bk{e)x^ + Bk+i{e)x^'^^ 


where Ao(e) is diagonal with distinct eigenvalues. A further normalization ofn—1 
non diagonal monomials in the numerator as in Lemma \3.1C\ can bring the system 
to a unique form. 
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Proof. We realize the modulus of such a family in the form (6.1) by Theorem 5.1 


Then the initial family and (6.1) are analytically equivalent since they have the 
same modulus. □ 

In the special case when the connection at e = 0 is irreducible, as noted above, 
we have given a new proof of the following theorem by Rostov: 

Theorem 6.3. f|12| ) A germ of family of linear differential systems unfolding an 
irreducible irregular non resonant singular point of Poincare rank k is analytically 


equivalent to a polynomial form as in (5.1) 

(6.2) Pe{x)y’ = (Ao(e) + Bi{e)x H-1- Bk{e)x'") ■ y, 

with Ao(e) diagonal with distinct eigenvalues. A further normalization ofn—1 non 
diagonal monomials as in Lemma \3.10\ can bring the system to a unique polynomial 
form. 


Proof. We realize the modulus of such a family in the form (6.2) by Theorem 5.1 


Then the initial family and (6.2) are analytically equivalent since they have the 
same modulus. □ 

Remark 6.4. The number of parameters in Theorem | Ai3| is optimal as remarked 
by Rostov in US]. Indeed, for each e, the modulus is described by 

• {k + l)n formal invariants (eigenvalues at each singular point) for each e; 

• 2k normalized upper or lower triangular Stokes matrices each with 


i(n—1) 


nontrivial upper or lower triangular entries. The equivalence defined above 
subtracts n — 1 parameters (the conjugation action of the scalar matrices 
being trivial) for a total of (n — l){kn — 1) parameters. 

All together, this yields kn'^ + 1 parameters. Now, each system (5.1) is described for 
each e by kn^ + n coefficients. This form is unique up to the action of a diagonal 
matrix, which allows further scaling of n — I coefficients. 


The number of parameters in Theorem \5.1\ can be explained in a similar way. 
The realized system p^(x)(x — R) = A„{x) depends on (fc + l)n^ + n parameters, 
which are the coefficients of A,:(x) (remember that ^^(O) is diagonal). A diagonal 
normalization reduces this number to (k + l)n^ + 1. On the other hand, 

• The Stokes matrices and the formal monodromy at the zeroes of p^ together 
give kn^ + 1 coefficients, quotienting by our equivalences, as above. 

• Moreover, the coefficients of the monodromy matrix around oo were 
chosen, basically arbitrarily in a small set. 

This leads to the same total of {k + l)n^ + 1. Hence, the number of parameters is 
explained by the full generality we introduced in the monodromy at infinity. 
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